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THE  STRUCTURE  OF  EDDIES  IN 
TURBULENT  FLAMES  -  I 

Department  of  Chemical  Engineering  and  Fuel  Technology 
University  of  Sheffield,  Sheffield  Si  3JD 

Norman  A.  Chigier  and  Andrew  J.  Yule 

Abstract 

This  report  includes  a  description  of  the  objectives  of  the 
investigation,  the  experimental  and  data  processing  techniques  and 

'a 

initial  results.  High  frequency  response  measurement  techniques 
are  used  with  data  processing  and  conditional  sampling  to  measure 
large  eddy  structure  and  coherence  in  the  transitional  flow  region 
of  turbulent  diffusion  flames.  Comparisons  are  made  between 
nonburning  and  burning  flow  conditions;  these  show  physical 
similarity  for  cold  and  burning  flows  in  the  transitional  region. 

Flow  visualization  shows  quasi- laminar  interface  layers  which  are 

distorted,  stretched  and  convoluted  by  vortex  interactions. 

^  0 ' 

Visualization  by^high-speed  Schlieren  colour  movie  show*  axisymmetric 
wave  and  vortex  ring  formation.  Velocity,  temperature  and 
ionization  density  measurements  show  that  the  transitional  flow  and 
the  potential  core  of  the  jet  are  much  longer  in  the  flame  than  in 
the  nonburning  jet.  A  marked  difference  is  found  between  the 
initial  fundamental  instability  frequencies  of  the  shear  layers; 

Tfie  most  amplified  frequency  predicted  by  stability  theory  for  cold 
inviscid  flow  is  foundin' the  cold  jet  but  not  the  flame.  Increase 
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in  viscosity  due  to  heat  release  causes  delay  and  damping  of  the 
processes  of  vortex  formation  and  coalescence;  dilatation  and 
buoyancy  forces  also  contribute  to  differences  between  the  flame 
and  cold  jet.  large  eddies  are  shown  to  play  an  important  role 

in  the  transitional  region  of  turbulent  jet  flames. 

J 
'  • 


Nomenclature 


Figure  1  shows  the  coordinate  system  in  the  experimental 
investigation  of  an  axisymmetric  diffusion  flame  in  a  uniform 
secondary  stream.  The  notation  below  covers  the  essential 
quantities  arising  in  data  acquisition  and  processing.  Further 
quantities,  occurring  for  example  in  various  methods  of  conditional 
sampling,  will  be  defined  as  they  arise. 

b(x)  mixing  layer  ’shear’  thickness,  b  =  -(U1-Ua)/Qu/3rjnav 

c_j  instantaneous  local  mass  fraction  of  specie  'j* 

0  diameter  of  burner  nozzle 

2 

Fr  Froude  number:  Fr  =  U_/gD 

J 

f  frequency  (Hz) 

fv  vortex  passing  frequency  in  transition  region 

I  intermittency  (with  value  1,  in  turbulent  fluid,  and  0,  in 

irrotational  fluid) 

I_._  instantaneous  saturation  ion  current 

SAT 

M  momentum  flux 

mass  flow  rate 

n  number  of  samples  in  one  time  history 

N  number  of  signal  time  histories 

P  instantanious  pressure 
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probability  density  function  of  a  (PDF) : 

p#a\  „  (Total  time  which  signal  is  between  g  and  da  ^ 

l  T  da  J [Limit  as] 

|da  -+■  0  J 

(p  can  be  obtained  by  ensemble  averaging  N  time  histories) 

cross  correlation  function  (autocorrelation  when  i  -  £) : 

T 

Rit(T)  =  limit^f  ai(t)ai(t+T)  dt 
T  JQ 

cross  correlation  coefficient: 


*U<T' 


Rit,Tl '  ‘'n*0'  v°» 


Reynolds  number :  Ret 
Reynolds  number :  Re 
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v/v 

V/V 


radial  distance  from  nozzle  centreline 

Strouhal  number  based  on  spectrum  peak  or  vortex  passing 
frequency:  StD  -  fpEAK  D/Uj  7  ^  0 


Strouhal  number:  St, 


fPEAK  5J^UJ 


Cross  (or  Auto)  power  spectral  density: 
-i  2irfx 


Sif<f) 


Sif(f) 


1  Ru  e_i 


^  x?  x„ 
T  i  «, 


i  i} 


dr;  and,  for  large  T: 

/  > 

/A 


length  of  signal  time  history 
time 

times  of  triggering  events  for  conditional  sampling 


instantaneous  velocity  components  in  x,  y  and  z  directions 
(W  *  0,  and  in  general  U  »  V) 

components  of  fluctuating  velocity: 
u  -  U-U;  v  =*  V-V  and  w  *  w 

(.Complex)  Fourier  transform  of  a^t) 

coordinates  based  on  nozzle  lip  (see  Fig.  1) 

coordinates  based  on  arbitrary  point  P,  e.g.  trigger  probe 
position 
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general  notation  for  a  dependent  scalar  or  vector  variable 


Y2(f) 


a  conserved  scalar,  e.g.  ft  is  the  local  concentration 
of  fluid  originating  from  the  burner  nozzle 

2  2 

coherence  function:  y  =  (S  )  /S 

U  ' 


ner  noz: 

,  *  J 


- 

nozzle  boundary  layer  shear  thickness:  J 

*j  ■  -oj  /<8  * 

x  =  0 

nozzle  boundary  layer,  displacement  and  momentum  thickness 

h  normalized  radial  coordinate  for  mixing  layer: 

n  =  y/x 


0 

v 

P 

T 


temperature  (K) 
kinematic  viscosity 
density 

probe  time  constant  or  time  relative  to  t  in  correlations  etc. 


Subscripts 


1  centre  line  value 

c  convection  velocity 

D  based  on  jet  diameter 

F  fundamental  frequency 

i,£,m  tensor  components 
I  irrotational  flow  value 

j  "jth"  species 

J  primary  (jet)  flow  value  (at  x  =  O) 

O  value  at  x  =  O 

S  secondary  flow  value 

T  turbulent  flow  value 

Superscripts 


denotes  a  fluctuating  component  measured  relative  to  a 
time  mean 

denotes  complex  conjugate 
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1 •  Introduction 

This  report  describes  the  experimental  and  data  analysis 
approaches,  both  in  use  and  under  development  by  the  authors,  to 
study  the  roles  of  'coherent  structures'  in  axisymroetric  turbulent 
gas  diffusion  flames.  The  experimental  results  from  the  first  17 
months'  work  on  the  research  project  are  described.  These  results 
concentrate  on  a  description  of  the  transitional  flame  structure 
near  the  burner  nozzle  and  comparisons  between  the  physical  structures 
of  cold  and  burning  transitional  flows. 

In  recent  years,  there  has  been  considerable  interest  in  the 
structures  of  large  eddies  in  many  types  of  turbulent  flow;  obser¬ 
vations  in  jet-type  free  shear  flows  are  of  most  interest  here.  The 
experiments  in  two-dimensional  mixing  layers  by  Roshko  ,  Browand  and 
their  coworkers  yielded  the  dramatic  results  shown  by  their  clear 
visualizations  of  repetitive,  coherent  eddies  which  appeared  to 
dominate  the  flows.  Observations  have  also  been  made  in  round  jets, 
for  example  by  Yule  and  Browand  and  Laufer  ,  which  have  indicated 
the  repetitiveness,  coherence  and  importance  of  the  large  eddies. 
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The  philosophy  behind  the  present  investigation  can  be  stated  in  its 
simplest  form:  Given  that  large  eddies  are  such  important  features 
of  nonreacting  shear  flows,  are  they  also  important  in  reacting 
flows  when  there  is  significant  heat  release  and,  if  so,  what  are 
the  structures  of  their  eddies  and  how  do  they  interact  with 
combustion? 

We  have  chosen  to  study  the  burning  round  jet  configuration, 
which  has  been  the  subject  of  many  previous  investigations  in  its 
cold,  nonreacting  form.  An  axisymmetric,  fuel  rich  gaseous  jet 
mixes  with  surrounding,  room  temperature,  air  and  burns  as,  essentially, 
a  diffusion  flame.  The  notation  used  is  shown  in  Fig.  1.  This 
burning  flow  will  be  referred  to  as  the  Jet  Flame,  and  the  nonburning 
flow  will  be  termed  the  Cold  Jet.  In  keeping  with  previous  inves¬ 
tigations  of  cold  jet  structure,  great  care  is  taken  to  ensure  that 
the  initial  conditions  at  the  nozzle  are  as  laminar  as  is  practically 
possible,  with  a  uniform  initial  velocity  profile.  Rather  surprisingly, 
there  has  been  little  attention  paid  in  the  past  to  the  effects  of 
initial  and  boundary  conditions  on  the  structure  of  turbulent  com¬ 
bustion  flow.  Reacting  turbulent  flow  is  often  implied  to  exist, 
per  se  ,  in  some  quasi-uni versal  fully-developed  turbulent  form, 
generally  referred  to  as  a  turbulent  diffusion  flame.  There  is  a 
striking  contrast  with  the  many  fundamental  investigations  of 
'classical'  cold  turbulent  shear  flows  in  which  initial  and  boundary 
conditions  are  controlled  with  great  care,  with  the  knowledge  that 
these  conditions  can  affect  the  turbulence,  and  its  development,  for 
some  considerable  distance  downstream.  Thus,  the  present  objective, 
to  obtain  a  fundamental  understanding  of  turbulent  flame  structure, 
with  emphasis  on  the  roles  of  the  large  eddies,  requires  that  the 
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jet  flame  Is  not  disturbed  significantly  by,  for  example,  pilot 
flames,  swirl,  free  stream  turbulence  or  recirculation  regions  used 
for  the  promotion  of  flame  stabilization.  In  order  to  stabilize 
the  flame,  low  velocities  are  required  and,  thus,  the  possibility 
arises  of  buoyancy  forces  becoming  significant  in  certain  regions 
of  the  flames.  However,  this  is  considered  to  be  a  small  price  to 
pay  for  the  achievement  of  simple  boundary  conditions  with  the 
resulting  simplification  of  analysis  and  data  interpretation  and 
the  production  of  a  'classical1  flow  which  can  be  readily  repeated 
and  reproduced. 

It  has  been  argued  that  the  study  of  such  a  well  controlled 
'simple'  flow  is  'pathological'  compared  to  most  practical  combustion 
flows  which  are  more  complex,  with  three-dimensionality,  confinement, 
swirl  and  recirculation.  The  authors  strongly  criticize  this  point 
of  view:  Until  we  have  achieved  a  fundamental  physical  understanding 
of  certain  basic  turbulent  reacting  shear  flows,  with  simple  boundary 
conditions,  we  cannot  hope  to  fully  understand,  and  reliably  model 
and  predict,  extremely  complex  situations.  However,  knowledge  of 
the  general  structure  of  such  complex  flows  is,  of  course,  important 
in  itself.  It  is  undeniable  that,  at  the  present  time,  there  is 
not  a  basic  understanding  of  even  'simple'  turbulent  flame  structure. 
This  is  obvious  from  the  multifarious  physical  models  and  analytical 
approaches  to  the  prediction  of  combustion  flows,  all  of  which  have 
one  common  factor:-  the  lack  of  sufficient  accurate,  time  resolved 
data  in  a  form  which  allows  conclusive  validation  of  the  assumed 
models  and  confirmation  of  modelling  predictions. 

It  is  argued  here  that  one  of  the  keys  to  the  understanding  of 
turbulent  combustion  lies  in  the  recognition  of  the  importance  of 
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strong  coherent  large  eddies .  Such  eddies  are  clearly  present  in 
many  cold  flows.  Furthermore,  examination  of  films  of  the  internal 
flows  in  practical  combustors,  large  and  small,  often  gives  indic¬ 
ations  of  coherent  burning  regions  moving  for  considerable  distances 
through  the  flow.  It  is,  thus,  reasonable  to  attempt  to  examine  a 
combustion  flow  from  a  quasi-La jrangian  viewpoint,  following  an 
eddy  from  its  formation  until  it  can  no  longer  be  recognized.  In 
the  process,  experimental  techniques  must  be  used  to  measure  how 
these  eddies  interact  with  each  other  and  how  they  interact  with 
the  processes  of  combustion.  For  example,  how  do  these  eddies 
influence  the  entrainment  and  mixing  of  gases  and  does  spatial  and 
temporal  coherence  in  such  eddies  result  in  'flamelets',  or  coherent 
regions  of  reaction,  which  are  recognizable  for  distances  comparable 
with  the  total  flame  length?  These  aims  require  measurement  and 

data  processing  techniques  of  the  multiprobe  conditional  sampling 

2  3 

type,  which  have  been  used  successfully  in  cold  flows  '  .  Accurate, 
multipoint,  time  resolved  measurements  in  flames  require  sophisticated 
diagnostic  techniques,  together  with  a  considerable  development 
effort,  as  is  described  in  Section  2  of  this  report. 

In  the  present  experiments  flame  structure  is  investigated  in  a 
logical  manner,  by  starting  at  the  burner  nozzle  and  extending  the 
detailed  measurements  gradually  downstream.  The  growth  of  large 
vortex- like  eddies  in  the  initial  transitional,  viscosity  dominated, 
mixing  layer- type  flow  is  first  examined.  The  interaction  and 
breakdown  of  these  initial  eddies  to  form  the  turbulent  flame,  and 
the  structures  and  convection  of  the  eddies  in  the  turbulent  flame, 
are  then  studied.  The  techniques  used  must  provide  quantitative 
data  on  the  structures  and  coherence  of  these  eddies,  and  also  on 
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the  interactions  between  combustion  and  the  eddies;  for  example, 
by  comparing  cold  and  burning  flows  of  the  same  jets,  and  by 
detecting  the  regions  of  reaction  associated  with  the  eddies. 

The  Regions  of  a  Jet 

It  is  useful  to  outline  current  knowledge  of  the  structure  of 
cold,  round  turbulent  jets  developing  from  laminar  nozzle  boundary 
layers.  As  sketched  in  Fig.  1,  the  flow  geometry  for  a  free  jet 
(i.e.  Ug  =  0)  can  be  divided  into  three  regions:  (i)  The  mixing 
layer,  or  potential  core,  region,  which  extends  up  to  x/D  =  5 
approximately.  This  is  characterized  by  a  central  region  of 
irrotational  flow  surrounded  by  the  high  shear,  mixing  layer; 

(ii)  A  development  region,  in  which  the  mixing  layer  turbulence 
interacts  and  adjusts  to  produce  the  third,  downstream  region; 

(iii)  the  fully  developed  region,  typically  beyond  x/D  =  15,  in 
which  all  mean  and  turbulence  quantities  have  approximate  (but 
never  exact)  self-preservation  so  that  the  mean  velocity  has  the 
similarity  form  0  =  x  (M  /p)  f  (— ) ,  where  x  is  measured  (in  this 

J  X 

region)  from  a  virtual  origin. 

In  addition,  there  is  a  very  important  division  between  the 
transition  region,  in  which  viscous  forces  are  influential,  and  the 
turbulent  region.  The  length  of  the  transition  region  is  dependent 
upon  the  initial  conditions  at  the  nozzle,  and  it  has  been  suggested 
that  xT  -  lOOO  6**,  although  this  has  not  been  shown  conclusively. 
Depending  upon  these  conditions,  the  transition  region  can  extend 
for  a  small  distance  along  the  mixing  layer,  xT/D  *0.25,  or  it  can 
extend  beyond  the  end  of  the  mixing  layer  region,  xT/D  >  5.  The 
initial  part  of  the  present  investigation,  which  is  described  in 
this  report,  lays  particular  emphasis  on  understanding  the  transitional 
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region  of  the  jet  flame,  and  making  comparisons  between  this  region 
for  cold  and  burning  jets. 

The  Importance  of  Transitional  Flow 

There  are  three  important  reasons  why  the  transitional  flow 
structure  is  of  importance  when  considering  the  structure  of 
diffusion  flames: 

(i)  Cold  flow  experiments  show  that  the  well  organized  vortical 

structures  in  transitional  flow  have  basic  similarities  with  the 

large  eddies  in  the  turbulent  region  of  the  same  flow,  further  down- 

3  4  5 

stream.  Indeed,  experiments  '  '  indicate  that  turbulent  eddies  in 
the  round  jet  evolve  directly  from  the  interactions,  coalescing  and 
breakdown  of  the  transitional  vortex  rings  formed  near  the  nozzle. 

(ii)  Viscous  forces  increase  by  an  order  of  magnitude  due  to  heat 
release,  so  that  transition  distances  in  flames  must  be  expected  to 
extend  much  further  than  the  distances  for  the  equivalent  nonburning 
flows. 

(iii) It  can  be  argued  that,  for  cold  nonreacting  flows,  the  details 
of  the  initial  conditions  and  the  transitional  flow  are  essentially 
forgotten’  beyond  some  distance  downstream,  beyond  which  Reynolds 
number  similarity  applies.  However,  this  argument  cannot  be  used 
for  burning  flows:  Heat  release,  mixing  and  product  and  pollutant 
formation  in  the  transition  region  are  integral  to  the  whole  flame 
and  are  not  ’forgotten’  by  the  turbulence  structure  after  the  end 
of  the  transition  region.  In  qualitative  terms,  it  is  useful  to 
think  of  the  transitional  flow  structure  as  being  determined  by 
initial  conditions  at  the  burner  nozzle,  while  the  turbulent  flow 
structure,  in  the  flame,  has  a  structure  determined  by  conditions 
existing  at  the  end  of  the  transition  region. 
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In  addition,  the  well  defined  growth  of  orderly  vortical 
structures  expected  in  the  transitional  jet  flame,  and  the  quasi- 
laminar  interface  burning  regions  which  should  be  associated  with 
them,  make  the  transitional  flow  an  ideal  environment  for  the 
application  and  development  of  the  conditional  sampling  and  flow 
visualization  techniques  used  in  this  investigation. 

The  'Physical  Structure  of  the  Round  Jet 

Yule^  used  flow  visualization,  time  and  frequency  domain 
processing  of  velocity  fluctuations  and  conditional  sampling  tech¬ 
niques  to  derive  quantitative  and  qualitative  knowledge  of  the 
mixing  layer  region  of  round  jets  at  different  Reynolds  numbers. 
Figure  2  sketches  the  physical  structure  which  was  found  at  all 
Reynolds  numbers  investigated  (9000,  <  ReD  ;2  x  lO5)  . 

Natural  instability  of  the  initial  laminar  shear  layer,  formed 
from  the  nozzle  boundary  layer,  produces  a  periodic  array  of  vortex¬ 
ring-like  vorticity  concentrations.  The  initial  stages  of  this 
instability  can  be  predicted  satisfactorily.  As  these  vortex  rings 
move  downstream,  they  generally  coalesce  with  neighbouring  rings,  so 
that  the  scale  and  separation  of  the  vortex  rings  increase  with 
distance  from  the  nozzle.  Although  point  velocity  measurements  in 
this  transitional  flc*?  exhibit  periodic  correlations  and  peaked 
spectra,  there  is  considerable  random  variation  in  the  strengths 
and  movements  of  the  coalescing  vortices.  Furthermore,  the  vortex 
rings  lose  their  phase  agreement  across  the  jet  as  they  move  down¬ 
stream.  Fluctuations  of  the  w  velocity  component  increase  with 
distance  from  the  nozzle  due  to  an  almost  linear  growth  of  orderly 
wave  deformations  of  the  cores  of  the  vortex  rings.  From  one  up 
to  three  regions  of  vortex  coalescence  have  been  observed  in  the 
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transition  region  depending  upon  the  nozzle  conditions,  which  also 
determine  vortex  passing  frequencies  and  length  of  the  transition 
region. 

Yule^  has  argued  that  the  last  coalescence  of  vortices,  prior 

to  turbulent  flow,  involves  rings  which  have  core  deformations 

larger  than  a  critical  size.  The  rings  entangle,  thus  producing 

enhanced  vorticity  stretching  and  small  scales  of  motion.  The 

remains  of  these  entangled  vortices  can  be  detected  in  the  turbulent 

region,  by  conditional  sampling,  up  to  5D  downstream  of  transition, 

and,  individually,  they  can  be  discerned  visually  for  several 

diameters  further  than  this,  i.e.  well  downstream  of  the  potential 

core.  These  modified  vortex  rings  thus  form  large  eddies  in  the 

turbulent  region  of  the  jet.  It  does  not  seem  to  be  valid,  in 

general,  to  consider  that  these  rings  completely  break  down  or 

disintegrate  into  turbulent  flow;  rather,  the  large  scale  vorticity 

concentrations  which  they  produce  retain  appreciable  coherence 

beyond  the  end  of  the  transition  region.  The  details  of  the  flows 

in  the  turbulent  mixing  layer  and  downstream  jet  regions  are  the 

subject  of  some  controversy.  There  appears  to  be  agreement  that 

the  large  eddies  are  strong,  coherent,  and  rather  similar  to  vortex 

rings  in  cross-section  when  measurements  are  ensemble  averaged  for 

3 

many  consecutive  eddies.  Yule  ,  however,  proposed  that  it  is  not 
reasonable  to  describe  these  eddies,  in  the  turbulent  region,  as 
vortex  rings  because  of  a  lack  of  circumferential  coherence  and  the 
smaller  scales  of  motion  within  them.  On  the  other  hand,  it  has  been 
proposed  that  the  eddies  of  the  turbulent  region  are,  in  fact, 
vortex  rings  but  their  circumferential  coherence  is  concealed  by 
'phase  scrambling'  effects  in  multipoint  measurements.  The  effects 
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of  phase  scrambling  in  the  present  investigation  are  discussed  in 
Section  3.  There  is  some  visualization  evidence,  but  no  quantit¬ 
ative  evidence,  that  the  turbulent  eddies  coalesce,  rather  than 
disintegrate  and  reform,  similar  to  the  eddy  coalescing  observed 
in  two-dimensional  mixing  layers  '  .  Furthermore,  the  basic 
large  eddy  structure  seems  to  be  similar  throughout  the  turbulent 
jet,  from  the  mixing  layer  and  into  the  fully  developed  region. 

These,  and  other  points,  await  elucidation  by  carefully  designed 
experiments . 

We  have  not  considered  here  the  evidence  of  large  eddy 
structure  in  two-dimensional  mixing  layers,  much  of  it  originating 
from  work  funded  by  Project  SQUID.  Although  one  can  expect 
general  similarities  between  the  two-dimensional  and  round  jet 
mixing  layers,  particularly  near  the  nozzle,  no  presumptions  are 
being  made  on  similarities  between  two-dimensional  and  axisymmetric 
flews.  Other  factors,  such  as  the  effect  of  a  secondary  stream 
on  the  jet  structure,  will  be  discussed  as  they  arise. 

Diffusion  Flams  Structure 

The  existence  of  vortex  rings,  near  the  burner  nozzle,  in 

flames,  has  been  recognized  for  many  years.  Consider,  for  instance, 

the  studies  of  Tyndall^  of  the  effect  of  sound  on  flames.  Grant^ 

has  described  some  observations  of  these  vortex  rings  and  Ballantyne 
0 

and  Bray  have  recognized  the  existence  of  such  vortex  structures  in 
certain  of  the  flames  which  they  have  studied.  However,  there  is  a 
lack  of  knowledge  of  the  growth  and  interactions  between  these  rings 
and  combustion,  and  even  such  basic  quantitative  information  as 
passing  frequencies,  coalescing  positions,  etc.  Essentially,  there 
is  no  systematic  information  on  the  interaction  between  combustion 
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and  transition  in  the  axi symmetric  flame.  For  the  downstream 

fully  developed  'turbulent'  flame  (the  assumption  of  local  turbulent 

flow  has  been  made  implicitly  in  the  past  without  proof)  there  are 

a  number  of  investigations  measuring  various  mean  and  rms  quantities. 

Often  these  have  indicated  the  importance  of  unmixedness,  which  is 

related  to  intermittency ,  probability  density  functions  and,  thus, 

9 

to  the  large  eddy  structure.  The  work  of  Kent  and  Bilger  and 
their  coworkers  is  perhaps  the  most  thorough,  and  the  most  relevant 
to  the  present  study.  However,  there  is  an  almost  complete  absence 

of  information  on  large  eddy  structure  in  turbulent  diffusion 
flames  in  the  published  literature  and  no  previous  investigations 
have  concentrated  on  the  roles  of  eddy  formation  and  coherence. 

2  Experimental  Facility  and  Techniques 

2 . 1  Introduction 

To  meet  the  objectives  described  above  requires  high  frequency 
measurement  techniques  with  good  spatial  resolution.  In  order 
that  conditional  sampling  and  cross-correlation  techniques  can  be 
applied,  to  educe  eddy  structures  and  their  convective  properties, 
facilities  must  be  developed  to  acquire  data  simultaneously  from 
several  probes  within  the  flame.  The  three  major  measurement  tech¬ 
niques  in  this  investigation  have  been  chosen  because  of  their 
relatively  advanced  stages  of  development,  and  relatively  good 
reliability,  resolution  and  accuracy.  The  design  development  and 
application  of  these  techniques  is  described  in  Sections  2.3,  2.4 
and  2.5  for  measurements  of  velocity,  temperature  and  ionization 
levels.  Ultimately,  to  quantify  the  flame  structure  fully,  a 
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technique  for  measurement  of  specie  fluctuations  is  required. 

Laser  Raman  techniques  may  eventually  fulfil  this  need,  but  versions 
of  more  conventional  species  concentration  measurement  techniques 
can  provide  useful  data,  as  will  be  described  in  Section  2.6.  It 
is  well  established,  from  cold  flow  studies,  that  flow  visualization 
plays  an  essential  part  in  any  attempts  to  study  the  coherence  and 
convective  properties  of  structures  in  quasi-Lagrangian  reference 
frames.  The  applications  of  flow  visualization  techniques  in  the 
present  investigation  is  described  in  Section  2.7. 

2.2  Facility  for  Production  of  Axisymmetric  Diffusion  Flame 

Figure  3  shows  a  schematic  diagram  of  the  apparatus  used  for 
the  production  of  an  axisymmetric  flame  in  a  coflowing  airstream. 

The  main  design  requirements  were  that  both  the  primary  jet  flow 
and  the  secondary  flow  had  low  turbulence  levels  with  uniform 
velocity  distributions .  To  enable  accurate  Laser  Doppler  Anemo¬ 
meter  (LDA)  measurements  of  velocity,  it  was  essential  that  the 
flame  was  surrounded  by  a  low  velocity  secondary  stream,  rather 
than  by  ambient  air.  This  is  because  both  the  primary  (fuel)  jet 
and  the  entrained  surrounding  air,  must  be  seeded  equally  with 
small  pair  tides  so  that  biasing  effects  are  minimized  in  the  velocity 
measurements.  All  experiments  reported  here  are  conducted  with  Ug  = 

0.73  m/s,  except  for  some  results  obtained  with  U  =0  (the  free  jet 

s 

or  free  jet  flame)  for  comparative  tests.  This  uniform  low  velocity 
secondary  flow  also  effectively  protected  the  flame  from  small 
draughts  in  the  laboratory,  thus  permitting  an  open  flow  with  easy 
access  by  probes.  The  secondary  flow  was  supplied  by  a  fan  via  a 
plenum  chamber  containing  sound-absorbing  foam,  to  ensure  that  fan 
and  pipe  noise  did  not  affect  the  flame.  The  primary  jet  is  fed 
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via  a  mixing  chamber,  diffuser,  settling  chamber  and  a  contoured 

nozzle  with  16:1  area  ratio.  The  mixing  chamber  has  three  inlets 

for  the  fuel  gas  (propane) ,  diluting  air  containing  LDA  seeding 

particles  and  nitrogen.  Nitrogen  was  used  only  for  safety,  to 

flush  the  system  before  and  after  runs;  however,  future  studies 

of  the  effects  of  Nj  dilution  are  envisaged.  All  flows  are 

measured  and  monitored  by  calibrated  rotameters.  The  facility 

permits  the  production  of  flames  ranging  from  'pure'  diffusion 

flames  (no  primary  air)  to  fully  premixed  flames .  Most  of  the 

initial  experiments  have  used  a  propane/air  primary  flow  with  an 

equivalence  ratio  of  10  (Jet  Flame  I) ,  which  is  sufficiently  fuel- 

rich  for  burning  to  occur  as  a  diffusion  flame.  This  propane /air 

mixture  has  been  used  to  achieve  a  reasonably  high  Reynolds  number, 
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ReD  =  1.04  x  10  ,  without  an  excessive  drain  from  the  propane 
bottles .  Volume  flow  rates  are  relatively  high  in  these  experi¬ 
ments  because  of  the  large  burner  nozzle  diameter  used;  D  =  25.4  mm. 
The  large  diameter  was  selected  to  enable  detailed,  accurate  measure¬ 
ments  of  initial  conditions  and  in  the  transitional  flow  near  the 
nozzle.  The  secondary  flow  convergence  is  square  in  cross-section 
and  has  an  exit  diameter  of  0.4  m.  Visible  flame  lengths  from 
1.5  m  to  2  m  are  found  and  the  initial  experiments  concentrate  on 
the  first  metre  of  flow.  The  complete  nozzle  and  secondary  flow 
assembly  can  be  traversed  radially  and  vertically  relative  to  the 
fixed  LDA  optics.  The  region  0  <  x<:  40  D  can  be  probed  by  the 
LDA.  Other,  intrusive,  probes  are  carried  by  traversing  apparatus 
mounted  on  vertical  lathe  beds  outside  the  secondary  flow. 

A  major  design  problem  with  the  apparatus  was  the  apparent 
incompatability  of  the  required  low  initial  turbulence  levels  with 


the  requirement  that  both  jet  and  secondary  flows  were  heavily 
seeded  for  the  LDA.  It  was  required  to  maximize  the  seeding  of 
the  flows  to  maximize  the  frequency  range  of  velocity  measurements. 
At  the  high  particle  densities  used  (of  the  order  of  1  particle 
per  mm3)  it  was  found  that  some  agglomeration  of  the  SiO^  particles 
was  inevitable.  Thus,  particles  of  the  order  of  2  ym  in  diameter 
were  detected,  although  the  nominal  particle  diameter  was  0.1  ym. 
Consequently,  any  fine  mesh  screens,  placed  in  either  the  jet  or 
secondary  flow  supply  systems,  rapidly  clogged  with  particles, 
causing  unacceptable  flow  variations.  After  several  different 
nozzle  configurations  were  tried,  the  problem  was  overcome  for  the 
jet  flow  by  omitting  all  smoothing  screens,  except  for  a  honeycomb, 
after  the  mixing  chamber  and  relying  on  long  lengths  of  low  velocity 
pipe  flow  to  permit  the  decay  of  the  turbulence  created  in  the 
mixing  chamber.  The  final  design  was  successful,  as  shown  by  the 
0.4%  turbulence  level  measured  at  the  nozzle  exit.  The  secondary 
flow  was  seeded  by  injection  of  particles  at  low  velocity  from 
holes  in  a  toroid  surrounding  the  nozzle  settling  chamber,  and 
downstream  of  the  screens  in  the  secondary  flow  settling  chamber. 

2 . 3  Laser  Doppler  Anemometer 

The  general  arrangement  of  the  LDA  system  is  shown  in  Fig.  4. 

A  one  Watt  Argon  Ion  Spectra  Physics  164  laser  is  used,  on  the  blue 
line,  A  =  0.488  ym,  with  a  TSI  beam  splitter  and  Bragg  cell  module. 
The  EMI  photomultiplier,  with  a  narrow  band  interference  filter, 
collects  light  in  the  forward  direction  and  slightly  off  axis.  The 
narrow  vertical  slit  in  front  of  the  photomultiplier  aperture  is  at 
the  focus  of  the  measurement  control  volume  by  the  collecting  lens. 
This  results  in  a  very  short  length  of  the  control  volume  being 


17 


I- 


4.7  mm 

2° 

6.1° 

599  mm 
5.87  pm 

lOOO  Hz 

±  1% 


viewed  by  the  photomultiplier.  Characteristics  of  the  LDA  system 
are: 

Beam  intersection  angle  =  4.8° 

2 

Measurement  control  volume  i/e  width  =  310  ym 

Length  of  measurement  control  volume 
seen  by  the  photomultiplier 
Typical  light  collection  angle 
Off-axis  angle  of  collection  optics 
Beam  focusing  less  focal  length 
Fringe  spacing  in  measurement  volume 
Typical  particle  arrival  rate 
(at  U  a  5  m/s) 

Accuracy  of  individual  velocity 
measurement  (U  “  5  m/s) 

The  LDA  operates  in  the  single  particle  counting  mode  so  that 
no  more  than  one  particle  at  a  time  is  permitted  in  the  region  of 
the  measurement  volume  viewed  by  the  photomultiplier.  A  typical 
signal  is  sketched  in  Fig.  5.  The  LDA  signal  processing  system 
has  been  designed  at  the  University  of  Sheffield  and  developments 
have  been  made  specifically  to  meet  the  requirements  of  the  present 
investigations.  Fringe  counting  logic  circuits  are  used  but 
various  interchangeable  circuit  boards  have  been  developed  to  meet 
different  requirements,  which  are  not  always  compatible,  such  as 
extremely  high  accuracy,  fast  data  rates,  or  the  uninterrupted 
acquirement  of  a  large  bulk  of  signals.  Figure  6  shows  the  logic 
block  diagram  for  the  processor  system  used  for  the  velocity 
measurements.  In  this  processor,  the  squaring  operations  on 
fringe  cycles  and  signal  envelope,  required  to  produce  gating 
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signals  for  the  operational  amplifiers ,  are  achieved  by  saturating 
amplifiers,  rather  than  by  Schmidt  triggers  as  have  been  used 
previously . 

The  counting  of  signal  fringe  cycles  commences  when  the 
signal  envelope  exceeds  a  threshold  level,  which  can  be  preset  to 
suit  the  prevailing  signal  quality  and  data  rate  requirements. 

The  processor  counts  6  cycles,  and  measures  the  count  duration  as 
the  number  of  cycles  of  a  20  mz  clock.  A  count  of  12  cycles 
is  also  made,  with  the  duration  T^  measured  by  a  10  MHz  clock. 
These  clock  frequencies  are  chosen  to  match  the  flow  velocities 
encountered;  higher  velocities  caul  be  measured  by  using  faster 
clocks.  The  D.E.C.  PDPll-03  minicomputer  acquires  T.  and  T,  _ 
for  each  signal  and  the  machine  code  acquisition  program  accepts 
the  data  if  T&  and  T  2  are  within  5%  of  each  other.  This 
comparison  criterion  value  was  derived  from  systematic  experiments 
which  showed  that  values  less  than  5%  did  not  significantly 
improve  data  accuracy .  There  are  two  main  sources  of  counting 
error.  The  first  is  due  to  the  counting  of  either  too  few  or  too 
many  fringe  cycles,  or  false  cycles,  due  to  noise  or  signal 
quality  variation  (e.g.  due  to  occasional  large  seeding  particles) . 
It  is  this  possible  error  which  is  recognized  and  thus  avoided  by 
making  a  comparison  test  between  the  times  for  two  numbers  of 
fringes.  The  second  general  source  of  error  arises  when  it  is 
not  possible  to  measure  exact  numbers  of  fringe  cycle  periods. 

There  are  several  possible  reasons  for  this;  for  example,  fringe 
cycles  do  not  all  have  the  same  amplitude  so  that  counts  can  begin 
and  end  at  different  positions  relative  to  cycles.  lhe  double 
fringe  count  comparison  can  remove  some  of  the  sources  of  error. 


19 


However,  this  type  of  error  is  minimized  only  by  maximizing  the 

fringe  quality  of  signals  and  carefully  tuning  the  electronics. 

It  was  necessary  to  acquire  only  the  two  clock  counts,  T  and 

6 

to  derive  a  particle  velocity.  A  system  developed  earlier, 
also  counted  the  actual  numbers  of  fringes  measured,  because 
signals  can  always  be  found  which  are  not  long  enough  to  satisfy 
the  specified  numbers.  However,  it  was  found  that  counts  of 
say  6  and  13  cycles,  rather  than  6  and  12,  are  recognized  by  the 
comparison  test  and  rejected.  Despite  the  precautions  taken, 
it  is  inevitable  that  small  amounts  of  poor  quality'  data  pass 
through  the  system  and  are  stored  by  the  computer,  together  with 
the  accurate  velocity  measurements.  These  'poor  quality'  measure 
ments  can  be  recognized  and  rejected  by  the  high  level  program 
used  to  further  process  the  date  for  determination  of  the 
velocity  statistics.  For  each  separate  measurement  location  in 
the  flow  a  file  of  velocity  measurements  is  first  presented  as  a 
probability  density  function  (PDF)  of  velocity.  Velocity  measure 
ments  which  are  inaccurate  enough  to  seriously  affect  the  accuracy 
of  velocity  statistics  are  recognized  by  their  occurrence  in 
extended  'tails'  at  both  ends  of  the  PDF.  Upper  and  lower  cut¬ 
offs  are  applied  to  reject  these  measurements,  which  account  for 
less  than  3%  of  the  total  data.  If  this  high  level  validation 

routine  is  not  followed,  errors  of  up  to  50%  can  occur  in  measure- 
— 2  *j  — 

ments  of  (u  )  ,  although  errors  in  U  are  insignificant.  These 
'tails'  in  LDA  velocity  PDFs  are  present  in  most  LDA  systems 
(including  commercial  systems)  known  to  the  authors.  In  general, 
their  users  are  aware  of  the  phenomenon,  but  it  is  worthwhile  to 
stress  that  large  errors  can  arise  if  due  care  is  not  taken  with 


data  validation. 

The  signal  processor  also  contains  a  'Real  Time'  Clock,  which 
times  the  occurrence  of  the  beginning  of  each  signal.  This  data 
is  also  acquired  by  the  computer  so  that  a  file  of  velocity  as 
a  function  of  time  is  stored.  Provided  data  rates  are  high, 
compared  with  the  most  energetic  frequencies  in  the  turbulence, 
accurate  time  averaged  quantities  can  thus  be  derived,  rather 
than  simple  ensemble  averages  of  signals. 

The  present  investigation  requires  a  maximum  error  of  no 
more  than  1%  in  individual  velocity  readings,  plus  high  frequency 
real  time  data  acquisition,  for  detailed  measurements  of  turbulence 
structure  and  to  enable  conditional  sampling  to  be  used.  The 
LDA  system  was  developed  and  checked  by  using  fibres,  attached  to 
rotating  disks,  traversed  at  known  velocities  through  the  measure¬ 
ment  volume,  and  also  by  directly  comparing  LDA  measurements 
with  hot  wire  measurements  in  the  same  flows. 

2.4.  Fine  Vive  Thermocouples 

The  types  of  thermocouple  used  in  this  investigation  are 
shown  in  Fig.  7.  Both  75  jim  diameter,  Pt-Pt/13%Rh  and  25  ym  diameter 
Pt/20%Rh-Pt/40%Rh,  probes  are  used.  The  latter  probes  are  fragile 
but  offer  a  higher  frequency  regime.  Unlike  resistance  probes 
which  suffer  from  insurmountable  drift  problems  in  flames,  the 
fine  wire  thermocouple  has  the  advantage  that  the  thermoelectric 
properties  of  the  junction  are  well  known  and  do  not  vary  with 
time.  When  relatively  thick  wires  are  used,  there  is  a  thermal 
lag  between  the  junction  temperature  and  the  local  gas  temperature 
in  a  fluctuating  temperature  field. 

To  the  first  order,  the  gas  and  wire  temperatures  are 
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usually  assumed  to  be  related  by  a  time  constant  so  that 


T,  ’  T„*T^  <« 
The  value  of  t  is  a  function  of  the  local  gas  properties  and 
velocity,  as  well  as  the  geometry  of  the  thermocouple.  The 
authors*-0  have  developed  a  method  of  measuring  x  with  the  thermo¬ 
couple  in  the  flame,  without  requiring  assumptions  to  be  made 
regarding  the  local  gas  characteristics  or  the  wire  geometry. 

As  shown  in  Fig.  8,  the  time  constant  measurement  system  is 
controlled  by  a  EEC  PDP  11-03  computer,  which  can  also  be  used 
subsequently  to  acquire  a  time  history  of  the  temperature  T^  which 
is  corrected  by  using  Eq.  (1) .  A  square  wave  overheating 
current  is  applied  to  the  thermocuple  and  a  portion  of  the  temper¬ 
ature  decay  curve  of  the  thermocouple  is  measured  after  each  over¬ 
heating  pulse.  As  shown  in  Fig.  9,  these  decay  curves  are 
ensemble  averaged  by  summing,  typically,  512  realizations.  By 
this  means,  flame  related  ten^serature  fluctuations  are  averaged 
out  to  give  a  smooth  decay  curve  which  describes  the  regime 
characteristics  of  the  thermocouple  (plus  its  electronics) .  The 
computer  program  matches  this  decay  curve  to  the  inverse 
exponential  decay  implied  by  Eq.  (1) ,  so  that  a  local  value  for 
x  is  derived.  In  addition*-0,  the  precise  shape  of  the  measured 
decay  curve  provides  information  on  the  deviation  of  the  regime 
characteristics  from  the  first  order  assumption.  In  fact,  the 
Fourier  transform  of  the  decay  curve  provides  the  true  response 
characteristics  of  the  probe.  It  was  found10  that  the  major 
source  of  deviations  from  the  first  order  equation  is  the  complex 
heat  transfer  characteristics  when  the  thermocouple  junction  bead 
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is  much  larger  than  the  wire  diameter.  For  the  25  vim  thermo¬ 
couples  used  in  this  study,  time  constants  of  lo  ms  are  common, 
and  accurate  signal  compensation  up  to  2  kHz  is  possible  by  using 
the  first  order  compensation  equation. 

Figure  10  shows  a  length  of  time  history  before  and  after 
compensation.  Compensation  (i.e.  signal  correction)  is  applied 
digitally  to  the  signal  by  high  level  program,  rather  than  by 
specially  designed  analog  circuitry  as  has  been  used  in  the  past. 
This  digital  compensation  procedure  is  convenient  and  flexible; 
the  main  requirement  is  that  there  is  very  little  electronic 
noise  introduced  in  the  signal  in  the  preamplifier  and  ADC 
electronics.  Such  noise  would  be  amplified  by  applying  Eq.  (1) 
to  the  digitized  signal.  Low  noise  is  achieved  by  the  use  of 
battery-driven  electronics  and  shielded  cables. 

2.5.  Ionization  Probes 

The  ionization,  or  Langmuir,  probes  used  in  this  investigation 
are  based  upon  the  designs  used  by  Lockwood  and  Odidi11.  As 
shown  in  Fig.  11,  the  probes  are  water-cooled  and  the  materials 
and  design  are  chosen  to  minimize  leakage  currents  and  flow 
interference.  In  hydrocarbon  flames,  the  ion  density  near  the 
reaction  zone  is  several  orders  of  magnitude  higher  than  else¬ 
where.  In  the  circuit,  shown  in  Fig.  12,  a  negative  voltage  is 
applied  to  the  platinum  electrode,  causing  a  current  to  flow  which 

is  dependent  upon  the  attraction  of  the  ions  to  the  electrode  (and 

12 

the  repulsion  of  electrons)  .  Several  theories  relate  the 
instantaneous  ion  current  to  the  ion  density,  flow  velocity,  probe 
geometry  and  biasing  voltage.  During  the  initial  stages  of  the 
present  investigation,  the  ionization  probe  is  being  used  to 
measure  the  flame  front  fluctuations .  Use  is  made  of  the  high 


frequency  response  (several  kHz) ,  good  spatial  resolution 

(<  0.5  mm^)  and  high  sensitivity  of  the  probe  to  the  local 

occurrence  of  reaction.  As  described  in  Section  3,  these 

characteristics  make  ionization  probes  ideal  for  measurement 

of  the  topography  of  reaction  zones  and  their  convective 

properties.  Further,  when  used  in  conjunction  with  probes 

measuring  other  quantities,  they  enable  conditional  sampling 

to  be  used  to  relate  averaged  local  eddy  structure  to  local 

burning  characteristics.  The  signals  are  amplified  by  low 

noise,  battery-driven,  preamplifiers  and  a  16  bit  ADC,  and  then 

fed  to  the  DEC  computer. 

2.6.  Other  Measurement  Techniques 

In  addition  to  the  main  techniques  described  above,  DISA 

hot  wires  and  constant  temperature  anemometers  have  been  used 

to  investigate  the  cold  jets.  Other  techniques  which  are 

available,  but  which  have  not  yet  been  applied  to  the  flames, 

include  streamlined  microphones,  for  the  measurement  of  near 

field  pressure  fluctuations,  and  standard  quartz  microprobe  (sonic) 

sampling  apparatus,  for  measurement  of  the  local  time  mean  species 

concentrations.  Small  modifications  to  the  LDA  system  will  permit 

the  measurement  of  conserved  scalar  fluctuations,  similar  to  the 

9 

methods  used  by  Bilger  .  Preliminary  LDA  experiments  have  also 
been  conducted  in  which  only  one  of  the  gas  streams  is  seeded. 

This  provides  ’weighted'  velocity  data  which  may  be  interpreted 
in  broadly  the  same  way  as  intermittency  conditioned  data  in  cold 
flows. 
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2.7.  Flow  Visualization 


Flow  visualization  is  considered  to  be  the  backbone  of 
this  investigation.  Direct  evidence  is  provided  of  'structure* 
in  the  flow  and  important  clues  are  provided  to  aid  in  the 
interpretation  of  quantitative  measurements.  The  main  flow 
visualization  system  is  shown  in  Fig.  13.  The  large  mirrors 
permit  a  20  D  length  of  flame  to  be  filmed  by  using  a  Hadland 
'Hyspeed'  camera,  which  has  a  maximum  framing  rate  of  20,000 
frames/second.  A  Schlieren  system  is  used,  as  interferometry 
would  provide  little  advantage  for  complex  three-dimensional 
turbulent  flows.  Colour  Schlieren  light  stops  are  used  to 
provide  a  clearer  qualitative  picture  of  the  different  density 
gradient  regimes.  Thus,  for  example ,  high  temperature  reacting 
interfaces  can  be  identified  by  streaks  of  a  certain  colour 
(typically  white) .  Kodak  Eastman  Ektachrome  7239  colour  cine 
film  is  used. 

No  attempt  has  been  made,  at  this  stage,  to  interpret  the 
films  in  terms  of  density,  or  temperature,  variations.  However, 
many  different  quantitative  analyses  can  be  applied  to  obtain 
important  information  on  the  development  and  interaction  of 
coherent  structures.  The  films  can  be  analysed:  (i)  manually, 
using  an  'editor';  (ii)  by  projection  on  a  screen  with  subsequent 
semi-automated  analysis  by  using  an  XY  coordinate  data-logger;  or 
(iii)  by  using  a  Quantimet  image  analysis  computer.  Films  are 
also  taken  directly  of  the  flame,  so  that  regions  of  luminous 
burning  are  seen  most  clearly.  The  differences  between  these 
films  and  the  Schlieren  films  are  striking,  as  will  be  discussed 
in  Section  4.  Signals,  from  probes  in  the  flame,  can  be 

25 


/ 


recorded  simultaneously  on  the  cine  films,  to  enable  comparisons 
to  be  made  between  observed  structures  and  the  signals  which 
they  produce. 

2.8.  Data  Acquisition  and  Processing  Facilities 

In  general,  all  data  is  acquired,  on-line,  by  DEC  PDP  11-03 
and  PDP-8E  computers  via  purpose  built  preamplifier,  ADC  and 
multiplexing  interfaces.  For  the  LDA,  for  example,  up  to 
40,000  signals  can  be  loaded  on  to  the  disk  of  the  PDP-8E  with 
the  PDP  11-03  acting  as  a  data  buffer  between  the  larger  computer 
and  the  LDA  processor /interface .  Fortran,  Basic  and  Machine 

Code  programming  are  used  for  data  processing.  16  bit  or  12  bit 
words  are  used,  depending  on  the  accuracy  required. 

A  spectroscope  Model  5D335  Real  Time  Analyser  is  used  to 
derive  power  spectra  from  continuous  signals. 

3  Data  Analysis  and  Interpretation 

3.1.  Laminar ,  Transitional,  Turbulent  and  Buoyancy  Influenced 
Regimes  of  Flames 

In  the  interpretation  of  quantitative  measurements  and  flow 
visualizations,  the  local  relative  importance  of  various  regimes 
of  flow  structure  must  be  recognized.  In  fact,  the  relative 
importance  of  these  various  regimes  can  only  be  deduced  from 
the  measurements,  and  flow  visualizations  themselves.  It  is 
important  to  avoid  an  assumption  of,  say,  local  turbulent  flow, 
without  supporting  data.  The  first  flames  to  be  studied  are 
basically  diffusion  flames  developing  from  relatively  turbulence- 
free  initial  nozzle  conditions.  It  is  found  that  ignition  and 
flame  occur  in  an  annulus  within  2  mm  of  the  nozzle  lip.  The 
initial  flame  is  almost  completely  steady  and  laminar,  both  in 
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this  region  and  also  for  several  mm  downstream.  However, 
although  the  reacting  region  of  the  flow  may  be  laminar,  it 
does  not  necessarily  follow  that  the  complete  width  of  the  jet 
is  laminar  in  the  initial  region.  Indeed,  it  has  been  found 
that  shear  layer  instabilities  are  already  developing  on  the 
inner,  high  velocity,  side  of  the  mixing  layer,  although  the 
outer  'sheath'  of  flame  is  completely  laminar . 

As  sketched  in  Fig.  14,  transitional  flow  develops  after 
the  laminar  initial  conditions.  By  definition,  in  the  tran¬ 
sitional  flow  viscous  and  inertial  stresses  are  of  the  same 
order  of  magnitude,  so  that  the  flow  is  Reynolds  number 

dependent.  Transitional  flow  is  complex,  even  in  the  cold  case; 

3 

however,  certain  repetitive  physical  stages  can  be  identified  . 

An  important  feature  is  the  development  of  three-dimensionality^ 
from  the  initial  axisymmetric  conditions,  although  the  dependence 
on  these  conditions  is  not  yet  fully  understood.  The  damping 
influence  of  viscosity  results  in  large  orderly  vortical 
structures,  with  little  evidence  of  smaller  eddy  scales.  The 
jet  flame  is  even  more  complex,  primarily  because  of: 

(i)  Increased  viscosity  due  to  heat  release.  The  relative 
importance  of  this  depends  upon  whether  this  viscosity  increase 
is  in  a  volume  of  fluid  where  viscous  forces  are  important. 

This,  in  turn,  depends  upon  a  complex  and,  as  yet,  unknown 
dependence  on  all  of  the  initial  and  boundary  conditions  of 
the  flame. 

(ii)  Dilatation,  or  volumetric  expansion,  which  must  also 
influence  the  transitional  flow  in  an  unknown  manner . 

The  initial  series  of  experiments  has  concentrated  on 
studying  the  interaction  between  combustion  and  the  fluid 
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mechanics  of  the  transitional  flow  with  the  objective  of  quan¬ 
tifying  these  effects. 

Further  downstream,  the  fully  turbulent  region  is  generally 
found,  although  for  low  velocity  flames  (or  premixed  flames) 
burning  can  be  completed  before  turbulence  has  developed.  Cold 
turbulent  flows  have  been  found^  to  have  instantaneous  three- 
dimensionality,  large  eddies,  a  wide  range  of  wave  numbers  in 
fluctuations  and  Reynolds  number  similarity .  It  remains  to  be 
shown  that  all  of  these  conditions  are  also  a  feature  of  the 
turbulent  flame.  As  with  the  initial  flow  near  the  nozzle,  it 
is  important  to  bear  in  mind  that,  while  most  of  the  central 
core  of  the  jet  flame  might  have  the  characteristics  of  turbulent 
flow,  this  is  not  necessarily  true  for  the  outer  low  velocity 
region  of  the  jet.  One  can  envisage  regimes  in  a  diffusion 
flame  where  heat  release  occurs  in  relatively  low  velocity,  outer 
regions  so  that  significant  reaction  can  occur  in  laminar  or 
transitional  flow  conditions,  while  the  bulk  of  the  jet  flame  is 
turbulent . 

The  occurrence  of  significant  buoyancy  forces  needs  to  be 
considered  when  interpreting  flame  data.  All  free  diffusion 
flames  contain  regions  where  forces  due  to  density  differences 
are  comparable  with  inertial  or  viscous  forces . 

Kennedy  and  Kent'*’’*  and  others  have  used  the  overall  Froude 

number  to  indicate  the  importance  of  buoyancy  forces.  They 

4 

showed  that  hydrogen  diffusion  flames,  with  FrD  =  7.5  x  10  , 
were  affected  by  buoyancy  forces  towards  the  end  of  the  flame 
where  'flapping1  was  observed.  The  various  independent  parameters 
which  determine  flame  structure  are  so  numerous  that  one  cannot 
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extend  this  observation  singly  to  the  present  experiments, 

which  use  a  larger  nozzle  and  propane  rather  than  hydrogen. 

Interest  is  being  initially  concentrated  on  the  first  20  D  of 

the  transitional  and  turbulent  flow,  where  buoyancy  forces  are 

of  much  less  relative  importance  than  at  the  end  of  the  visible 

flame,  at  x  =  80  D.  This  flame.  Jet  Flame  I,  has 
4  4 

ReD  =  1.05  x  10  ,  and  Fr^  =  1.63  x  lO  .  The  principal  para¬ 
meters  for  the  jet  flame  are  given  in  Table  1.  The  effects  of 
buoyancy  on  specific  regions  of  the  flame  will  be  examined  by 
calculating  the  relative  order  of  magnitude  of  inertial,  viscous 
and  buoyancy  forces  from  measured  data  of  velocity  and  temperature 
Buoyancy  forces  will  dominate  in  regions  of  high  temperature  and 
low  velocity.  Such  conditions  are  found  at  the  outer  edge  of 
the  jet  flame. 

3.2.  Space,  Time  and  Frequency  Domain  Analysis  of  Data 

An  ultimate  objective  of  this  investigation  is  to  measure 
coherent  structures  in  flames  in  a  quasi-Lagrangian  reference 
frame,  so  that  structures,  and  their  interactions,  are  measured 
as  they  convect  downstream.  This  objective  requires  specialized 
conditional  analysis  techniques  of  both  multipoint  data  and  flow 
visualization  films.  Conventional  time  average  measurements  also 
form  an  essential  part  of  the  investigation.  The  detailed 
information  on  coherent  structures  obtained  by  direct  measurement 
requires  to  be  related  to  structures  determined  by  classical 
statistical  averaging  techniques.  The  various  quantities  to  be 
measured  in  the  flames  (and  cold  jets)  are  listed  and  defined  in 
the  Notation.  The  single  point  time  histories  of  velocity, 
temperature  and  ionization  density  are  processed  digitally  for 


selected  regions  of  the  flows.  Power  spectra  and  autocorrelations 
are  derived  using  fast  Fourier  transform  techniques.  In 
addition,  more  detailed  mapping  of  each  flow,  including  the 
initial  conditions,  is  made  in  which  point  measurements  are 
processed  to  give  time  mean  and  rms  fluctuating  components  of 
the  various  quantities,  together  with  probability  density  functions. 
The  data  derived  from  the  simultaneous  multiprobe  measurements 
(currently  under  development  for  conditional  sampling) ,  will  also 
be  processed,  in  selected  regions  of  the  flame,  to  give  cross¬ 
correlations  or  cross-power  spectra.  Clasically  averaged  single 
and  multipoint  data  will  be  used,  as  far  as  possible,  to  derive 
values  of  individual  terms  in  the  partial  differential  equations 
for  continuity,  momentum  and  energy.  For  example,  simultaneously 
acquired  LDA  and  thermocouple  data,  measured  with  the  thermo¬ 
couple  close  to  the  LDA  measurement  volume,  will  enable  a  good 
approximation  to  Favre  averaged  velocity  data,  u2.  Ionization 
probe  and  thermocouple  data  will  be  processed  to  provide  infor¬ 
mation  on  local  reaction  rates. 

3.3.  Conditional  Sampling  for  the  Study  of  Flame  Structure 
Conditional  sampling  techniques  are  being  developed  for 
the  study  of  combustion  flows,  using  single  and  multipoint 
measurements.  The  expression  ’Conditional  Sailing’  is  used  in 
a  geieral  sense  to  distinguish  from  the  classical  statistical 
averaging  procedures  noted  in  3.2,  those  techniques  of  data 
analysis  which  use  ensemble  averaging  to  examine  flow  features  of 
interest,  with  the  averaging  conditioned  either  by  a  feature  of 
the  signal  itself  or  by  an  event  or  quality  of  another  signal 
acquired  simultaneously  elsewhere  in  the  flow.  A  wide  range  of 
conditional  sampling  techniques  has  been  applied  to  cold  turbulent 


flows.  Some  of  these  techniques  have  been  outlined  by  Davies 
and  Yule^,  Van  Atta^  and  Laufer^.  Yule1^  has  attempted  to 
bring  some  unification  to  the  conditional  sampling/coherent 
structures  field  by  including  a  systematic  description  of  the 
interpretation  of  various  conditioned  data  analysis  techniques. 

There  are  two  principal  'coherent  structures'  of  interest 
for  the  present  investigation: 

(i)  Interface  burning  regions  The  structures  and  coherence  of 

such  regions,  where  they  exist,  are  of  extreme  importance  from 

the  point  of  view  of  overall  flame  structure,  pollutant  emissions 

and  mixedness.  Physical  models,  such  as  the  'wrinkled'  flame 
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front  and  Spalding's  ESCIMO  concept,  have  been  used  in  the 
past  to  describe  premixed  flames. 

(ii)  Large  Eddies  including  the  'birth',  structure  and  inter¬ 
actions  of  these  eddies  and,  in  particular,  their  interactions 
with  the  interface  burning  regions. 

Many  cold- flow  conditional  sampling  techniques  are  not 
useable  in  reacting  flows  until  further  advances  are  made  in 
flame  measurement  techniques.  For  example,  LDA  techniques 
cannot  currently  be  used  to  detect  the  highest  frequencies 
occurring  in  the  turbulence  so  that  velocity  signals  cannot  be 
reliably  processed  to  obtain  intermittency  information.  On 
the  other  hand,  the  many  additional  dependent  variables  occuring 
in  reacting  flows  compared  with  cold  flows,  vastly  increases  the 
permutations  of  conditioning  criteria  and  conditioned  signals, 
which  could  be  used.  It  is,  thus,  necessary  to  be  selective, 
and  concentrate  on  designing  techniques  to  provide  direct  infor¬ 
mation  on  large  eddy  structures  and  the  interface  burning  region. 
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Some  of  these  procedures  are  noted  here: 


(i)  Conditioned  Single  Point  Data  The  ionization  probe  provides, 
perhaps,  the  most  useful  signal  for  the  application  of  ensemble 
averaging  techniques.  The  sharp  'spikes'  typically  found  in 
ionization  density  signals  are  produced  by  reaction  zones  and 
provide  a  clear  and  distinct  signal  of  reaction  occurring  at 
interface  layers  moving  through  the  flow.  By  fixing  triggering 
levels  in  signal  time  histories,  information  can  be  derived  on 
the  frequency  of  occurrence  of  these  layers,  on  the  thickness  of 
the  layers  and  on  variations  found  in  the  structures  of  reacting 
interface  layers  at  a  fixed  position. 

The  ten£>erature  and  velocity  time  histories  from  single  point 
measurements  are  examined  by  using  'pattern  recognition’  techniques. 
The  objective  of  using  these  techniques  is  to  determine  the 
proportion  of  a  signal  that  can  be  attributed  to  structures  which 
have  a  high  degree  of  local  repetitiveness.  Such  single  point 
data  cannot,  however,  provide  information  on  the  coherence  of 
such  structures.  Examination  of  time  histories  in  this  manner 
is  an  essential  prerequisite  of  multiprobe  conditional  sampling 
experiments,  as  information  can  be  deduced  concerning  the  triggering 
criteria  which  are  likely  to  be  most  usefully  applied.  Yule"* 
proposed  that  the  simultaneous  acquisition  of  single  point  signals, 
together  with  flow  visualization  films,  provides  the  optimum 
method  of  designing  conditional  sampling  techniques. 

(ii)  Multiprobe  Conditional  Sampling  of  Local  Structures  Signals 
acquired  simultaneously  by  two,  or  more,  probes  with  variable 
spatial  separation,  are  processable  to  obtain  averaged  local 
structures  of  the  eddies  of  interest.  Figure  15  shows  an  example 
of  a  situation  in  which  a  fixed  ionization  probe  is  used  to  detect 
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the  passing  of  a  reacting  interface  layer.  The  signals  from 
the  LDA  or  thermocouple  are  ensemble  averaged,  at  different 
positions,  when  the  ionization  signal  peaks  occur,  so  that  the 
average  temperature  and  velocity  fields  are  measured  relative 
to  the  interfaces.  Variations  and  elaborations  of  this  con¬ 
figuration  are  envisaged;  for  example,  two  adjacent  ionization 
probes  can  be  used  to  trigger  the  sampling,  with  the  'condition' 
that  passing  interfaces  are  at  a  specific  angle  to  the  axis. 
Signals  from  the  LDA  or  thermocouple  can  also  be  used  for  the 
triggering  criteria.  In  cold  jets3,  conditional  sampling, 
relative  to  potential  core  velocity  peaks,  reveals  averaged 
large  eddy  structures.  Figure  15  shows  a  situation  where  the 
burning  interface  encircles  the  eddy,  as  is  found  to  be  the  case, 
for  the  first  10  D,  of  Flame  I.  The  eddies  have  vortex-ring 
shapes  and  the  flow  is  similar  to  the  transitional  flow  in  cold 
round  jets. 

A  related  technique  utilizes  an  array,  or  rake,  of  probes 
stretched  radially  across  the  flow.  A  rake  of  8  thermocouples, 
spaced  at  10  sm  intervals,  has  been  constructed  for  this  purpose. 
This  arrangement  allows  conditional  averaging  of  the  radial 
temperature  distribution,  and  also  provides  a  means  for  searching 
for  repetitive  patterns  in  the  instantaneous  radial  tenperature 
distribution. 

(iii)  Measurement  of  Structure  Coherence  and  Convection  by 
Multiprobe  Conditional  Sampling  The  techniques  described 
in  the  previous  section,  for  measuring  local  structure,  can  be 
extended  to  measure  the  changes,  coherence  and  convection  down¬ 
stream  of  the  structures.  The  signal  recovery  probes  are 
positioned  at  different  positions  downstream  of  the  triggering 
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probe,  and  the  signals  are  averaged  with  a  time  delay  sufficient 
to  permit  the  convection  of  the  sampled  structures  to  the 
recovery  position. 

Figure  16  shows  a  situation  in  which  two  ionization  probes 
are  used  to  measure  the  convection  and  changes  in  dimensions  of 
reacting  interface  layers.  Results  obtained  from  these  measure¬ 
ments  will  be  compared  with  the  equivalent  space-time  cross¬ 
correlations  and  indications  will  be  sought  of  the  extent  to 
which  these  correlations  can  conceal,  or  ’smear’,  the  true 
coherence  of  convected  structures.  Extraction  of  information 
concerning  the  true  coherence  and  structure  of  an  event  from  the 
signal  properties  introduced  by  'phase  scrambling'  effects  are 
discussed  in  Section  3.5.  When,  for  example,  coalescence  of 
the  structures  of  interest  is  known  to  occur,  the  triggering  and 
conditioning  procedures  require  to  be  specifically  designed  to 
educe  this  information  from  the  recovered  signals . 

3.4.  Flow  Visualization:  Quantitative  Applications 

The  colour  Schlieren  and  direct  photography  high-speed  cine 
films  of  the  flames  are  processed  to  provide  quantitative  data 
on  flame  characteristics.  Such  quantitative  analyses  are  possible 
for  at  least  the  first  20  D  of  the  flow  in  'Jet  Flame  I',  where 
repetitive  structures  and  interfaces  are  very  clearly  visible. 

It  is  not  yet  known  how  far  this  quantitative  analysis  can  be 
extended  downstream.  Vortex,  or  instability  wave,  passing 
frequencies  can  be  measured  from  films  in  the  transitional  flow 
near  the  nozzle.  Both  in  this  region  and  further  downstream, 
the  properties  of  very  clear  interface  burning  regions  are  derived 
by  measuring  the  changes  in  shapes  and  separations  of  these 
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interfaces  as  they  are  convected  downstream.  Analysis  of  films 
provides  probability  distributions  of  the  shapes  of  interfaces 
at  different  positions  and,  as  shown  in  Section  5,  this  inform¬ 
ation  provides  data  on,  and  great  insight  into,  the  smearing,  or 
phase  scrambling,  inherent  in  multiprobe  point  measurements. 

3.5.  Data  Interpretation  in  Term  of  the  Coherence  of  Structures 
and  Phase  Soranbling  Effects 

In  transitional  and  turbulent  flows,  there  is  evidence  that 
structures,  either  eddies  or  interfaces,  retain  a  high  degree  of 
spatial  coherence,  or  identity,  as  they  are  convected  downstream. 
The  major  problem  in  space-time  cross-correlations  and  in  con¬ 
ditional  sampling  experiments  is  to  recognize  this  coherence . 

This  coherence  must  be  recognized  and  established  in  individual 
structures  undergoing  some  degree  of  change  in  dimensions, 
geometry  and  trajectory  as  they  move  downstream.  Other  possible 
phenomena,  such  as  three-dimensionality  in  structures  and  their 
movements  and  the  coalescing  of  structures,  combine  to  produce 
low  levels  of  recovered  signal,  or  cross-correlation,  in  spite 
of  the  coherence  of  individual  structures.  This  'smearing'  of 
multipoint  data  will  be  referred  to  as  'phase  scrambling',  since 
the  major  contributor  to  the  loss  of  recovered  signals  is 
produced  by  variations  of  the  phases  of  individual  sampled 
signals  resulting  from  combinations  of  the  various  phenomena. 

Figure  17  illustrates  the  phenomenon  for  the  case  when  a 
characteristic  signal  a^tx)  is  measured  at  times  relative  to  the 
passage  of  the  structure  past  the  trigger  probe  (signal  a^(t)) 
further  upstream.  For  the  flame,  the  signal  can  represent 
signals  from  ionization  probes  or  thermocouples  when  reacting 
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interface  layers  are  present.  The  conditionally  sampled  ensemble 
average  <  u2(t)  >  has  a  low  amplitude  compared  with  the  amplitudes 
of  individual  signals,  so  that  am  average  convection  velocity 
<Uc>  is  the  only  useful  quamtity  that  can  be  extracted  from  the 
ensemble  average.  Various  alternative  more  complex  methods  of 
ensemble  averaging  can  be  employed.  Figure  17  illustrates  one 
such  procedure,  where  the  phase  scrambling  effects  are  ’unscrambled1 
by  ensemble  averaging  the  recovered  signal  at  times  relative  to 
the  signal  peaks,  which  occur  at  for  the  nth  sampled  structure. 
This  procedure  provides  a  more  representative  ensemble  averaged 
structure,  with  a  high  amplitude  and  width  more  truly  represent¬ 
ative  of  the  individual  structures.  Additional  information  such 
as,  for  example,  the  probability  distribution  of  the  arrival 
times,  is  also  provided.  Such  a  procedure  is  only  possible 
for  signals  with  sharp  peaks  which  are  clearly  recognizable. 

Thus,  this  approximates  to  the  convection  of  one-dimensional 
structures.  More  complex  situations  require  more  complex  con¬ 
ditioning  techniques.  For  example,  rather  than  using  signal 
peaks,  the  phase  differences  between  individual  recovered  signals 
can  be  derived  by  phase  shifting  each  individual  signal  to 
obtain  the  best  fit  with  the  predetermined  ensemble  average 
<  >  .  A  new  'unscrambled'  ensemble  average  can  then  be 

derived  by  taking  these  phase  shifts  into  account. 

Particular  problems  are  encountered  when  velocity  signals 
are  used  for  conditional  sampling.  If  an  eddy  retains  coherence 
with  distance  downstream,  this  should  be  most  clearly  recognizable 
in  scalar  fields  associated  with  the  eddy  or  in  the  vorticity 
distribution  in  the  eddy.  The  local  fluctuations,  in  say  the 


— 


0  velocity  component,  will  not  necessarily  be  as  recognizably 
repetitive  with  increasing  distance  downstream  because  of  their 
interdependence  with  the  vorticity  and  potential  flow  fields 
surrounding  them. 

Figure  18  illustrates  the  more  general  situation  regarding 
coherent  structures.  Here  the  nth  structure  is  identifiable 
(indeed  it  is  defined  by)  a  three-dimensional  distribution  of 
some  scalar  or  vector  quantity  a(x"  ,  y"  ,  z"  )  .  The  tguasi- 
Lagrangian'  coordinate  system  is  based  upon  a  reference 
position  embedded  in  the  structure.  This  coordinate  system 
need  not  necessarily  be  identified  with  the  same  ' fluid  element' 
throughout  the  structure's  lifetime;  it  is  more  likely  to  be 
either  the  equivalent  of  a  centroid  of  the  structure  or,  alter¬ 
natively,  a  position  in  the  structure  which  produces  the  point 
signal  which  is  recognized  by  the  signal  conditioning  technique. 

In  the  same  way,  there  is  no  reason  why  a  coherent  structure 
should  be  identified  with  the  same  packet  of  fluid  throughout 
its  existence.  In  fact,  it  can  be  argued  that  this  is  manifestly 
not  the  case  for  most  coherent  structures.  Thus,  the  larger 
scales  in  the  early  transitional  regions  of  shear  flows  can  be 
considered  as  waves  moving  through  the  fluid  and  the  large  eddies 
in  the  turbulent  flow  continually  entrain  new  fluid,  and  also 


shed  old  fluid 


a  flow  of  fluid  across  the  interface 


interface  layers  require  a  continual  supply  of  fresh  mixture  for 


diffuse  into  the  irrotational  flow 


As  shown  in  Fig.  18,  the  structure  moves  downstream  until 


at  some  time  t,  it  can  be  identified  in  a  certain  position  and 

with  a  certain  change  in  the  shape  and  distribution,'  with 

a  =  a(x”  ,  y"  ,  z"  ,  t)  .  There  may  also  be  some  loss  of 

coherence.  Such  a  loss  may  involve  a  breakdown  in  the  structure 

so  that  it  is  no  longer  identifiable  as  such.  A  rather 

different  type  of  breakdown  occurs  if  the  structure  amalgamates 

with  a  neighbouring  structure.  The  precise  definition  and 

concept  which  describe  the  coherence  of  individual  eddies  are 

somewhat  arbitrary  and  even  subjective.  This  is  opposed  to 

the  coherence  functions  in  statistical  time  averaged  analyses, 

although  these,  and  the  equivalent  cross-correlations,  can  give 

a  very  poor  indication  of  the  true  physical  coherence  of 

individual  structures.  An  important  objective  of  current 

coherent  structure  studies  is  the  identification,  specification 

and  definition  of  'coherence'  as  applied  to  individual  structures. 

A  step  towards  this  understanding  will  be  made  when  the  phase 

scrambling  effects  can  be  quantified.  If  the  structures,  shown 

in  Fig.  18,  are  ensemble  averaged  over  many  realizations,  at 

time  t,  an  ensemble  averaged  structure  is  obtained.  <  a(x',  y',  z',x)>, 

where  the  coordinate  system,  (x',  y',  z’),  is  an  Euler ian  system 

based  upon  P,  the  most  probable  position  of  the  reference  point 

Pn  averaged  over  many  structures  at  time  t  =  x^.  The  ensemble 

averaged  structure  will  be  smeared  because  of  the  variations  is 

individual  structures.  Thus,  the  nth  structure,  shown  in  Fig.  18, 

is  at  a  position  (x^  ,  Yp  «  zp  )  relative  to  P.  There  will, 
n  n  n 

thus,  be  a  probability  density  function  describing  the  probability 
of  the  reference  point  of  a  structure  being  at  a  position 

(x',  y',  z')  at  t.  which  can  be  written  P„_,™,mr, (xj, ,  y',  z’  t,). 

P  P  P  1  STRUCTURE  P  JP  P  1 
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If  all  structures  at  time  x^  are  identical  so  that  signals 
at  t  =  t1  are  scrambled  only  due  to  shifting  in  the  positions  of 
structures  from  one  to  the  next,  phase  scrambling  is  a  function 
of  PSTRUCT  only#  s°  that  the  ensemble  averaged  structure  at 
t  =  t  is  related  to  the  a  field  of  an  individual  structure  by: 


<  a(x' ,  y' ,  z'  ,-r^  > 


PSTRUC  (XP'  yPf 
TURE 


ZP'  Tl> 


(2) 


a(x”  =  x*-x^,  y"  =  y’-yp,  z"  =  z’-z^,  x^  dx^  dy^  dz^ 


where  the  integral  is  taken  over  the  total  volume  in  which  the 
reference  point,  Pn»  of  structures  can  be  found  at  the  time  delay 
x  =  x^ 

Further  contributions  to  the  loss  of  signal  can  be  expected 

due  to  a  lack  of  repetition  in  the  triggering  criterion  relative 

to  the  structures  and  a  lack  of  repetitiveness  in  the  individual 

structures.  Further  forms  of  Equation  (2)  can  be  formulated  to 
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describe  these  effects.  It  is  then  necessary  to  design 
experimental  systems  to  quantify  these  contributions  and,  thus, 
attempt  to  extract  information  on  the  true  coherence  and  structures 
of  individual  eddies  from  the  data.  These  concepts  are  equally 
applicable  to  the  quantitative  analysis  of  flow  visualizations 
and  data  on  phase  scrambling,  derived  from  Schlieren  films,  as 
discussed  in  Section  5. 


4  Preliminary  Results 

The  initial  experiments  have  concentrated  on  the  first  20  D 
of  Jet  Flame  I  and  the  equivalent  nonburning  flow,  Cold  Jet  I. 

The  independent  variables  for  these  flows  are  specified  in  Table  1. 


4.1.  Flow  Visualization  of  Jet  Flame  I 

Figure  19  shows  still  photographs  of  the  complete  flame. 

From  visual  observations,  the  flame  appears  to  be  ’turbulent' 
beyond  x  =  5  D.  High-speed  photography,  with  close-up  views, 
reveals,  however,  that  the  basic  flow  (at  least  in  the  first 
2o  D  which  have  been  analysed  so  far)  consists  of  distinct  flow 
regimes  which  contain  coherent  structures.  A  qualitative 
description  of  this  flow  is  given  here,  while  Section  5  discusses 
quantitative  data  derived  from  the  films. 

Figure  20  is  a  single  frame  from  a  high-speed  film  of  the 
flame,  showing  the  region  8  x  4  16  D.  Figure  21  gives  a 
sequence  of  four  frames,  at  8  ms  intervals,  showing  the  region 
16  D  <  x  ^  24  D.  Both  of  these  films  were  obtained  by  directly 
filming  the  flame  so  that  light  areas  represent  regions  of  high 
flame  luminosity.  Figure  20  shows  the  appearance  of  azimuthally 
distributed  streaks  of  luminous  burning.  These  streaks  are 
themselves  distributed  fairly  periodically  with  longitudinal 
distance  along  the  jet  flame.  These  streaks  are  considered  to 
be  formed  by  the  azimuthal  wave  instability  of  the  vortex  rings 
in  the  first  part  of  the  jet  flame.  It  is,  thus,  deduced  that 
similar  three-dimensional  instability  mechanisms  can  be  found  in 
both  the  jet  flame  and  in  the  cold  round  jet.  The  overall 
structure  of  a  cold  jet  is  sketched  in  Fig.  2.  In  spite  of  the 
general  similarity,  quantitative  differences  must  be  expected 
between  the  cold  and  burning  cases,  and  this  will  be  shown  to 
be  the  case. 

The  sequence  of  frames  in  Fig.  21  shows,  on  the  left  hand 
side  of  the  jet  flame,  bulges  of  luminous  burning  which  appear 
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to  combine  and  produce  an  outward  moving  tongue  of  flame .  This 

is  evidence  of  the  coalescing  of  structures  in  the  jet  flame. 

This  phenomenon  could  be  seen  often  in  this  region  but  the  event 
is  far  from  repetitive  when  examined  using  flame  luminosity. 

It  is  found  that  flow  visualization,  by  simple  use  of  flame 
luminosity,  is  not  the  most  useful  technique  for  deriving 
qualitative  and  quantitative  information  on  flame  structure. 

There  are  two  main  reasons  for  this:  (i)  structures  in  three- 
dimensional  flows  are  hidden  by  luminous  areas  on  the  outside  of 
the  flow;  (ii)  the  occurrence  of  flame  luminosity  is  sensitive 
to  fluid  mechanical  and  chemical  kinetic  criteria.  Such 
criteria  may  only  be  met  in  certain  regions  of  structures,  so 
that  a  flow  visualized  by  flame  luminosity  may  appear  relatively 
disorderly  while,  fluid  mechanically,  quite  orderly  structures 
are  present.  The  use  of  colour  Schlieren  flow  visualization 
with  high-speed  cine  films,  indicated  a  substantial  orderliness, 
in  the  first  20  D  of  flow  at  least.  Figure  22,  showing  the 
region 0<x<10D using  colour  Schlieren,  indicates  this  orderliness. 
The  reaction  zone,  up  to  approximately  x  =  1  D,  appears  to  be 
completely  laminar  and  almost  cylindrical,  kut  with  a  larger 
diameter  than  that  of  the  nozzle.  Further  downstream,  beyond 
x  =  2  D,  periodic  waves  develop  in  the  reacting  interface.  These 

waves  are  believed  to  be  caused  by  the  entrainment  action  of  the 
developing  ring  vortices.  The  cores  of  these  vortices  are  found 
on  the  ’high'  velocity  part  of  the  mixing  layer  (see  Section  5) . 

The  high  temperature  interfaces  can  be  seen  to  stretch  around 
each  vortex;  discrete,  individual  interfaces  break  away  and 
separate  from  the  initially  continuous  interface.  Beyond  x  *  6  D, 
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an  array  of  'S'  shaped  hot  reacting  can  be  seen  in  the  central 
part  of  the  mixing  layer.  These  interfaces  and  their  associated 
vortices  or  eddies  are  seen  to  coalesce  in  the  region  4  D  <  x  <  8  D. 
Similar  'S'  shaped  interfaces  persist  beyond  the  region  of 
coalescence,  in  the  high-shear  region  of  the  jet  flame;  they 
have  been  found  at  least  up  to  x  =  20  D,  which  is  the  furthest 
distance  downstream  which  has  been  examined  so  far.  Beyond  the 
vortex  coalescing  region  (x  >  8  D) ,  the  shape,  changes  in  shape 
and  movements  of  the  interfaces  are  visually  less-repetitive 
than  further  upstream.  Some  individual  interfaces  can  be  seen 
to  undergo  extinction  of  reaction,  while  others  coalesce  with 
neighbours,  but  with  less  organized  repetition  than  in  the  region 
x  <  8  D. 

The  cine  films  show  two  fairly  distinct  regimes  co-existing 
in  the  first  20  D  of  the  flame.  There  is  the  inner,  fast  moving 
flow,  which  is  in  the  central  region  of  the  jet  flame, 

0  <  r  <  0.4  D,  and  an  outer  slower  moving  flow  which  contains 
significant  reacting  regions.  Observations^'4  in  cold  jets  have 
shown  that  the  interacting  vortices  and  eddies  in  the  mixing  layer 
shed  fluid  outwards  to  form  a  relatively  disorganized  sheath  of 
slow-moving  fluid.  It  seems  likely  that  this  same  mechanism  is 
also  present  in  the  jet  flame  and  is  responsible  for  the  outer, 
slower  moving  component  of  the  flow. 

All  the  phenomena  which  have  previously  been  observed  in 
the  mixing  layers  of  cold  jets  have  now  also  been  identified  in 
the  jet  flame.  Furthermore,  these  phenomena  have  distinct  and 
clear  effects  on  the  dynamics  of  coherent  interface  burning 
regions  which  appear  to  be  the  dominating  components  of  the 
reacting  flow.  Reaction  also  affects  the  shape  and  dynamics  of 
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the  vortices  and  eddies,  through  dilatation  and  increased  viscosity 
effects.  In  the  cine  films,  vortex-like  eddies  were  indentif iable 
up  to  at  least  x  =  20  D,  with  considerably  more  clarity  them  for 
cold  jets.  It  will  be  shown  that  there  are  important  differences 
in  the  rates  of  spread  of  the  mixing  layer,  the  vortex  frequencies 
and  the  transition  distance  between  the  jet  flame  and  the  cold 
jet. 

4.2.  Cold  Jet  I  Investigation 

Measurements  have  been  made  in  Cold  Jet  I  with  a  hot  wire 
anemometer  and  laser  anemometer  (LDA)  to  establish  very  definit¬ 
ively  the  initial  conditions  for  the  mixing  layer  in  the  plane  of 
the  jet  exit.  It  has  been  assumed,  but  will  subsequently  be 
verified  by  LDA  measurements,  that  the  initial  conditions  for  Jet 
Flame  I  are  the  same  as  for  Cold  Jet  I.  Velocity  measurements 
have  also  been  made  within  the  cold  mixing  layer  and  jet  so  that 
data  can  be  compared  directly  with  LDA  data  from  Jet  Flame  I. 

This  comparison  provides  information  on  the  interaction  between 
combustion  and  fluid  dynamics  in  the  transitional  and  turbulent 
flow. 

Figure  23  shows  velocity  profiles  measured  near  the  inner 
nozzle  lip  at  the  nozzle  exit  plane  (x/D  =  0.04)  for  the  flow 

conditions  of  Cold  Jet  I  and  also  at  a  lower  Reynolds  number. 

T  ~2  *j 

Measurements  of  (u  )  show  very  low  values,  with  (u  )  /U  no  more 

J 

than  0.005.  Thus,  the  initial  jet  flow  was  essentially  laminar. 
The  laminar  inner  boundary  layer  was  found  to  be  thicker  than  in 
most  jets  that  have  previously  been  investigated.  Figure  24 
shows  that  data  for  the  ’maximum  shear'  thickness  of  the  boundary 
layer  fit  the  equation 
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(3) 


S/D  =  9.5  Re  ~’3 

u  D 

In  previous  studies,  a  typical  value  for  the  proportionality 
constant  in  the  equation  is  3.  The  explanation  for  this 
relatively  thick  boundary  layer  is  found  in  the  special  require¬ 
ments  in  the  design  of  the  unusually  long  nozzle  and  entry  pipe 
to  provide  low  turbulence  levels  without  the  insertion  of  screens 
or  gauzes  and,  at  the  same  time,  to  prevent  separation  of  the  low 
velocity  secondary  flow  on  the  outside  of  the  nozzle.  It  can  be 
argued  that  the  inner  nozzle  boundary  layer  may  not  be  the  only 
factor  determining  the  initial  jet  instability  when  there  is  a 
secondary  stream.  The  boundary  layer  on  the  outside  of  the 
nozzle,  and  also  the  wake  caused  by  the  nozzle  lip,  could  con¬ 
tribute  towards  the  instability,  but  these  effects  should  be 

small  when  U  «  U  ,  as  is  the  case  for  Cold  Jet  I.  In  an 
S  J 

idealized  situation,  if  the  initial  velocity  deficit  at  the  lip 

is  quickly  filled  by  entrainment  from  the  secondary  stream,  an 

initial  shear  layer  thickness  b  can  be  defined  such  that 

o 

b  =  (1  -  U_/U_)  6,.  The  expected  variation  of  b  with  Re  . 

O  S  vJ  vJ  O  D 

when  U  is  kept  constant,  is  included  in  Fig.  24. 

b 

Figure  25  shows  radial  traverses  at  x  =  2  D  in  Cold  Jet  I, 
with  (U  /U  =  0.11)  and  without  (U„/UT  =  O)  the  secondary  flow. 

It  is  seen  that  the  velocity  gradient  is  the  central  part  of  the 
mixing  layer  is  unaffected  by  the  presence  of  the  secondary  flow. 
Small  differences  are,  however,  found  between  the  two  cases,  in 
the  potential  core  flow  of  the  jet.  It  is  thought  that  these 
small  differences  can  be  attributed  to  changes  in  the  vortex 
dynamics  in  the  initial  transitional  flow.  Figure  25  also  shows 
the  good  agreement  found  between  the  LDA  and  hot  wire  data.  The 
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differences  between  the  cases  U  /U  =  O  and  0.11  are  also  seen 

S  J 

— *r  h 

in  the  centreline  values  of  U  and  (u  )  shown  in  Fig.  26.  The 
rate  of  increase  of  fluctuations  in  the  potential  core  is  sig¬ 
nificantly  lower  for  the  case  U  /U,  =  0.11.  This  damping 
effect  may  be  considered  as  being  caused  by  a  'stretching'  of 
the  transitional  flow  by  the  secondary  stream.  Measurements 
show,  however,  that  the  secondary  stream  does  not  affect  the 
initial  instability  frequency;  nor  does  it  affect  the  frequency 
spectra  peaks  throughout  the  transitional  jet.  A  full  descrip¬ 
tion  of  the  effect  of  a  secondary  stream  on  the  transitional 
flow  in  round  jets  requires  further  study. 

At  this  stage  of  the  investigation,  evidence  suggests  that 
the  presence  of  the  low  velocity  secondary  stream  does  not 
drastically  affect  the  overall  physical  structure  of  the  jet; 
a  number  of  recognizable,  small  effects  have  been  detected.  All 
subsequent  data  refer  to  jets  and  jet  flames,  with  the  secondary 
flow. 

Comprehensive  mapping  of  amplitude  spectra  of  u  in  Cold  Jet  I 
has  been  made  to  derive  detailed  knowledge  of  the  transitional 
flow.  Figures  27  and  28  show  spectra  at  the  centreline  for 
0  <  x/D  <  6.  The  fundamental  instability  frequency  occurs  at 
215  Hz;  the  amplitude  increases  to  reach  a  maximum  at  x/D  =  3.5. 
This  'growth'  corresponds  to  the  rolling  up  of  the  initial  shear 
layer  instability  to  form  near -periodic  vortex  rings. 

The  subharmonic  of  the  fundamental  frequency,  peaking  at 
104  Hz,  is  just  discernable  at  x/D  =  1;  it  is  caused  by  the 
pairing  or  coalescing  of  vortices  as  they  move  downstream.  After 
x/D  =  3.5,  the  fundamental  frequency  peak  decays  until,  beyond 
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x/D  =5,  it  ts  lost  in  the  typical  turbulent  spectrum.  However, 
the  subharmonic  peak  increases  to  a  maximum  value  at  x/D  =  4  and 
it  is  still  discernable  at  x/D  =6. 

The  spectra  at  r/D  =  0.5,  shown  in  Figs.  29  and  30  are 
strikingly  different  from  those  on  the  centreline  (r/D  =  O) . 

The  fundamental  said  subharmonic  frequencies  are  the  same  as  those 
found  in  the  potential  core,  although  there  is  a  small,  but 
noticeable,  decrease  in  the  subharmonic  peak  frequency,  with 
increasing  x.  In  the  transitional  mixing  layer,  the  fundamental 
frequency  peak  is  greater  in  amplitude  than  the  subharmonic 
amplitude  until  both  peaks  disappear  at  x/D  =  5.  Between 
x/D  =  1.5  and  x/D  =  3,  there  is  a  very  rapid  increase  in  the 
amplitude  of  the  subharmonic;  this  coincides  with  the  region  of 
coalescing  of  initial  vortices.  The  mixing  layer  spectra,  in 
Figs.  29  and  30,  also  display  first  and  higher  order  harmonics 
of  the  fundamental  frequency.  These  higher  harmonics  are 
associated  with  the  viscous  core  of  the  vortices. 

Figures  31  and  32  show  spectra  at  different  radial  positions 
at  x/D  =  4  and  x/D  =  6.  The  peaks  that  are  evident  across  the 
complete  mixing  layer  at  x/D  =  4,  disappear  at  x/D  =  6,  apart  from 
the  subharmonic  peak  which  is  discernable  at  r/D  =0.  On  the 
basi..  of  the  establishment  of  turbulent  forms  for  spectra  inside 
the  mixing  layer,  the  transition  length  of  Cold  Jet  I  is  deduced 
to  be  xT  =  5  D.  By  coincidence,  this  roughly  corresponds  with 
the  end  of  the  potential  core. 

4.3.  Measurements  in  Jet  Flame  I 

Comprehensive  measurements  of  velocity,  temperature, 
ionization  current  and  other  quantities  in  flames  are  continuing 
to  be  made  and  periodic  reports  of  this  information  will  be  made. 


In  this  report,  a  selection  of  data  is  presented  with  the  aim  of 
providing  clear  definitions  and  descriptions  of  the  initial  and 
boundary  conditions.  The  data  and  flow  visualization  from  the 
flame  and  cold  jet  presented  in  this  report  also  provide  insight 
into  the  combustion-flow  interaction  in  the  initial  part  of  the 
flame. 

Figure  33  shows  LDA  velocity  profiles  in  the  initial  part 
of  Jet  Flame  I.  In  general  terms,  the  flow  in  the  flame  is 
similar  to  that  in  the  cold  jet.  There  is  a  central,  uniform 
velocity,  potential  core  with  an  outer  'cylindrical'  mixing 
layer.  The  mean  temperature  distribution  at  x/D  =  2,  which  is 
also  included  in  Fig.  33,  has  a  peak  at  the  outer  part  of  the 
mixing  layer.  The  location  of  this  temperature  peak  within  the 
mixing  layer  is  determined  by  the  location  of  the  stoichiometric 
fuel/air  mixture  ratio.  In  Fig.  34,  the  velocity  distributions 
at  x/D  =  4  are  compared  for  cold  and  burning  jets.  The  potential 
core  has  ended  for  the  cold  jet  while,  for  the  burning  jet,  the 
core  still  extends  out  to  r/D  =0.4.  The  rate  of  spread  of  the 
mixing  layer  into  the  potential  core  of  the  flame  is  very  much 
slower  than  in  the  cold  jet.  These  differences  can  be  attributed 
primarily  to  the  increased  viscosity  and  dilatation  effects  on 
the  initial  transitional  vortex  flow  in  the  flame. 

Flow  visualization  indicated  that  the  flame  was  stabilized 
in  a  steady,  laminar,  cylindrical  interface  which  was  attached  to 
the  nozzle  lip.  This  laminar  flame  appeared  to  be  steady  up  to 
x/D  =  2  and  had  a  diameter  greater  than  the  nozzle  diameter. 

Figure  35  shows  temperature  and  velocity  profiles  at  x  =  2  mm 
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in  the  flame  stabilization  region.  It  is  seen  that  the  peak 
temperature,  corresponding  the  the  centre  of  this  initial  laminar 
diffusion  flame,  occurs  0.2  D  outside  the  high  shear  region  of 
the  velocity  distribution.  The  temperature  increase  in  this 
high  shear  region  is  no  more  than  300  K,  but  the  viscosity 
changes  due  to  this  increase  might  be  expected  to  be  sufficient 
to  affect  the  initial  shear  layer  instability  growth.  The 
low  velocity  region  between  0.5  r/D  and  0.65  r/D,  where  flame 
stabilization  actually  occurs,  will  be  investigated  in  detail 
using  a  newly  acquired  Bragg  cell  frequency  shift  unit  for  the 
LDA. 

The  spread  of  the  mixing  layer  into  the  potential  core 

region  of  the  flame  is  greatly  reduced  compared  to  the  cold  jet. 

However,  on  the  outer,  low  velocity,  section  of  the  mixing  layer, 

where  high  temperatures  are  found,  velocities  are  higher  than 

those  found  in  the  cold  jet  mixing  layer.  These  increased 

velocities  are  primarily  due  to  dilatation,  although  buoyancy 

forces  will  also  provide  contributions  to  the  acceleration  of 

the  outer  flow.  For  example,  at  x/D  =  2,  in  the  high  temperature 

region  at  r/D  =  0.6,  the  velocity  in  the  flame  is  U/U  =  0.38, 

J 

while,  at  the  same  position,  in  the  cold  jet,  U/U  =  0.2.  The 

J 

T  h 

rate  of  increase  of  the  peak  'turbulence'  intensity  (u  )  is 
much  less  in  the  jet  flame  than  in  the  cold  jet  (Fig.  34) . 

In  the  first  10  D  of  flow,  the  ionization  probe  was  found 
to  give  significant  signals,  only  on  the  inner  side  of  the  mixing 
layer  of  the  jet  flame,  i.e.  for  r/D  <  0.5.  Figure  36  shows 
a  typical  spectrum  of  the  ionization  current  signal  at  x/D  =  8; 
r/D  =  0.5 .  Various  distinct  peaks  can  be  seen  in  the  spectrum. 
Some  of  these  peaks  can  be  identified  with  frequencies  derived 
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by  analysing  high-speed  cine  films .  It  can  be  seen  from 

Fig.  36,  that  the  fundamental  instability  frequency  of  the  cold 

jet  cannot  be  identified  in  the  spectrum  for  the  flame. 

5.  Discussion  and  Concluding  Remarks 

The  initial  conditions  of  the  cold  and  burning  jets, 

including  the  initial  shear  layer  instability  and  flame  stabil- 

19 

ization,  will  first  be  discussed.  Michalke's  theory  predicts 
the  most  amplified  shear  layer  instability  frequency  as  a  function 
of  the  initial  shear  layer  thickness.  Figure  37  shows  the 
predictions  from  this  theory  in  the  form  of  the  Strouhal  number 
St^  .  It  is  seen  that  the  measured  fundamental  frequency  of 
Cold  Jet  I  is  in  good  agreement  with  the  prediction  of  Michalke. 

It  is,  thus,  reasonable  to  assume  that  both  the  cold  and  burning 
jets  become  unstable  naturally,  so  that  there  is  no  influence 
from  the  low  initial  turbulence  levels  nor  from  external  sound 
sources.  Most  of  the  initial  turbulence  at  the  nozzle,  where 

T  ij 

(u  )  /U  <  0.4%,  is  concentrated  near  the  30  Hz  frequency,  as 
J 

can  be  seen  in  Fig.  27.  The  small  fluctuations  are  thought  to 

be  the  result  of  an  'organ- type'  resonance  of  the  nozzle  and  inlet 

pipe.  This  resonance  does  not  affect  the  natural  shear  layer 

instability  which  occurs  at  215  Hz. 

The  transition  distance  of  the  cold  jet,  deduced  from  the 

spectra  measurements  described  in  Section  4.2,  is  included  in 

Fig.  38.  Figure  38  also  includes  measurements  of  the  transition 

distance,  as  a  function  of  Reynolds  number  Re^  made  by  Yule^  in 

a  50.8  mm  diameter  cold  round  jet.  For  the  latter  jet,  the 

nozzle  boundary  layer  thickness  obeyed  the  proportionality 
-I. 

6j/D  =  3  ReQ  .  Thus,  for  a  given  value  of  ReQ,  the  50.8  mm 
diameter  jet  had  a  much  thinner  nozzle  boundary  layer  than  the 
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present  25.4  mm  diameter  jet.  It  is  seen  from  Fig.  38  that 

the  data  do  not  satisfy  the  proportionality  x  'v  6^  which 

20 

has  been  proposed  by  several  authors  .  In  addition,  no 
reasonable  agreement  can  be  found  with  the  alternative  suggestion 
that  the  transition  Reynolds  number,  x  U  /v,  is  approximately 

T  vJ 

constant.  It  is  clear  that  round  jet  transition  is  still  far 
from  being  understood  and  further  work  is  required  to  elucidate 
the  relative  importance  of  viscous  damping  (including  the  role 
of  the  initial  nozzle  boundary  layer) ,  coalescence  and  the 
growth  of  three-dimensionality .  It  seems  reasonable  to  assume 
that,  when  a  transitional  mixing  layer  flow  extends  towards  the 
end  of  the  potential  core,  interaction  across  the  jet  will 
enhance  the  growth  of  three-dimensionality,  and  thus  the  develop¬ 
ment  of  turbulent  flow.  Thus,  for  example,  if  ReD  is  fixed  and 
6^  is  gradually  increased  by  some  means,  xt  may  indeed  initially 

increase  almost  linearly  with  6,.  However,  when  x  becomes 

J  T 

comparable  with  the  potential  core  length,  transition  is 
accelerated  and  xt  then  increases  less  quickly  with  increasing 


A  further  important  consideration  is  the  role  of  vortex 
coalescence  in  transition.  A  feature  of  transitional  jet  flow 
is  the  near-periodic  vortex  ring  formation  and  coalescence  which 
occurs  in  well-defined  regions.  These  vortex  rings,  and  their 
coalescence,  are  recognized  by  single,  double  or  triple  spectra 
peaks.  The  present  Cold  Jet  I  has  only  one  recognizable 
vortex  coalescence  prior  to  turbulent  flow.  The  large  eddies, 
after  transition,  appear  to  coalesce,  but  they  lack  the  repetit- 
ativeness  required  for  the  production  of  spectra  peaks  and  sub¬ 
harmonics  .  Up  to  three  transitional  coalescence  regions  have 
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been  observed  in  various  jets  '  .  It  appears  that  the  latter 
stage  of  transition  involves  the  coalescence  of  vortex  rings 
which  have  developed  core  waves  greater  than  some  critical  size. 

It  can  be  surmised  that,  if  the  end  of  transition  is  identified 
with  a  particular  coalescence,  a  sufficiently  detailed  investig¬ 
ation  of  as  a  function  of  Re^  for  a  given  jet,  might  reveal 
’ripples’  in  the  distribution  corresponding  to  transition 
occurring  after  one,  two,  three  or  perhaps  more  coalescings. 

Figure  39  includes  the  axial  distributions  of  Strouhal 
numbers  corresponding  to  spectra  peaks,  and  thus  vortex  passing 
frequencies,  in  the  cold  jet.  Data  have  only  been  included  for 
regions  where  peaks  can  be  identified  in  the  spectra.  For 
the  jet  flame,  spectra,  such  as  those  shown  in  Fig.  36,  cannot 
yet  be  unambiguously  interpreted  in  terms  of  vortex  passing 
frequencies.  However,  in  the  first  20  D  of  the  jet  flame,  at 
least,  the  vortices,  eddies  and  interfaces  associated  with  them 
can  be  identified  sufficiently  clearly  in  Schlieren  cine  films 
for  passing  frequencies  to  be  measured  accurately.  Figure  39 
includes  Strouhal  number  distributions  derived  from  eddy  passing 
frequencies  measured  from  the  Schlieren  film.  The  differences 
between  the  cold  and  burning  jets  are  large  and  their  nature  was 
rather  unexpected.  The  initial  frequency,  identified  at  x/D  =  2, 
was  more  them  twice  the  initial  frequency  found  in  the  cold  jet. 
The  coalescence  of  these  initial  vortices  in  the  flame  occurred 
in  the  region  4  D  <  x  <  8  D  as  opposed  to  2  D  <  x  <  4  D  in  the 
cold  jet.  This  stretching ’  of  the  initial  part  of  the  flame 
transition  can  be  attributed  to  both  the  effect  of  dilatation  in 
accelerating  the  vortices  and  also  the  temperature-related 
viscosity  increases  which  decrease  the  rate  of  growth  of  the 
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vortices . 


Measurement  of  100  consecutive  vortices  between  3  D  <  x  <  5  D 

yielded  an  average  vortex  convection  velocity,  U  =  0.92  U  in 

c  J 

this  region  of  the  flame.  This  contrasts  with  U  =  0.6  U,  for 

c  J 

the  cold  jet.  Individual  vortices  in  the  flame  had  convection 
velocities  in  the  range  0.86  U,  <  U  <  0.98  U  .  Further 

downstream,  the  vortices,  or  eddies,  were  found  to  accelerate 
so  that  between  12  D  <  x  <  16  D  the  average  convection  velocity 
is  U  =  1.06  u  .  In  this  same  region,  individual  convection 

C  J 

velocities  were  measured  in  the  range  0.94  0  <U  <  1.22  U  . 

J  c  J 

Using  Schlieren  and  direct  visualization  techniques,  the 
appearance  of  the  initial  part  of  the  jet  flame  is  shown  in 
Fig.  40.  The  visualization  by  flame  luminosity  indicates 
predominantly  fuel-rich  burning  regions,  which  are,  initially, 
associated  with  specific  three-dimensional  features  of  the 
flame.  In  the  first  two  diameters  of  flow,  the  cylindrical 
laminar  flame  is  separated  from  the  inner,  high  shear  flow,  as 
is  shown  by  the  Lift  and  thermocouple  data  in  Fig.  35.  The 
velocity  profile  at  x  =  2  mm  in  the  high  shear  region,  resulting 
from  the  detached  nozzle  boundary  layer,  differed  insignificantly 
from  that  measured  in  the  initial  cold  flow  by  the  hot  wire 
anemometer.  It  cannot,  therefore,  be  argued  that  an  altered 
initial  velocity  profile  is  responsible  for  the  higher  initial 
instability  frequency  in  the  burning  jet.  Other  factors  that 
could  influence  the  initial  frequencies  are,  feed-back  pressure 
or  sound  fields  set  up  downstream  of  the  nozzle  by  combustion/flow 
interaction.  In  addition,  the  high  temperature  at  the  outside  of 
the  initial  shear  layer  could  cause  damping  of  instabilities, 


resulting  in  a  reduced  'effective  thickness*  of  the  shear  layer  compared 
to  that  found  in  the  cold  jet.  Further,  more  detailed,  measure¬ 
ments  require  to  be  made  in  the  initial  flow  region  in  order 
to  explain  fully  these  differences.  Further  insight  is  being 
sought  by  carrying  out  a  mathematical  analysis  using  a  finite 
difference,  time-stepping  model  to  investigate  the  effects  of 
combustion  on  shear  layer  instability  and  vortex  roll-up. 

Between  2  D  <  x  <  4  D  Schlieren  films  show  the  outer 
reacting  interface  being  entrained  by  the  inner  periodic  vortices. 
Significant  flame  luminosity  is  only  found  after  x  =  6  D  in 
the  region  of  coalescence  of  the  initial  vortices .  This 
luminosity  occurs  as  axial  streaks  around  pairs  of  coalescent 
vortices.  These  streaks  correspond  to  the  development  of  three- 
dimensionality  in  the  vortices,  similar  to  that  previously 
observed  in  cold  round  jets'?  After  x  =  8  D,  continuous  inter¬ 
face  burning  regions  begin  to  break  up  into  detached  'S'  shaped 
interfaces  associated  with  the  leading  edge  of  each  eddy.  In 
this  region,  although  there  is  some  noticeable  circumferential 
symmetry  across  the  flame,  it  is  debatable  whether  the  observed 
structures  can  be  labelled  'vortex  rings'.  Beyond  x  =  12  D 
the  movements  of  these  reacting  interfaces  become  noticably  less 
repetitive  although,  individually,  they  generally  remain  coherent 
up  to  at  least  x  =  20  D. 

A  striking  feature  of  the  jet  flame  downstream  of  approximately 
x  =  10  D  was  the  appearance  of  an  outer,  slow  moving  component  of 
the  flow.  Tongues  of  reacting  or  hot  gas  detached  themselves 
from  the  main,  inner  flow,  containing  the  eddies  and  interfaces. 

This  outer  flow  moved  downstream  at  a  velocity  of  approximately 
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0.3  U  ,  while  the  inner  eddies  moved  with  U  *  U  .  These 

J  c  J 

outer  tongues  of  flow  have  a  similar  appearance  to  the  three- 
dimensional  'jets'  of  relatively  disorderly  fluid  which  have 
been  observed  to  be  ejected  from  cold  jets. 

The  experiments  have  shown  that  reaction,  in  the  first 
20  D  of  the  jet  flame,  at  least,  is  concentrated  largely  in 
interfaces.  An  individual  region  of  reacting  interface  can 
generally  be  followed  from  its  formation  near  the  nozzle, 
through  coalescence,  up  to  x  =  20  D.  These  interfaces,  therefore, 
can  reasonable  be  termed  'coherent  structures'  in  the  flame,  and 
they  are  also  associated  with  coherent  vortex-like  eddies. 

One  objective  of  this  study  is  to  measure  the  structures  and 
histories  of  these  coherent  structures.  As  a  first  approach 
to  this  end,  a  conditional  sampling  technique  has  been  applied 
to  the  analysis  of  high-speed  cine  films.  The  shapes  of  fifty 
interfaces  were  recorded  at  the  times  when  the  passed  through  a 
'trigger'  position,  x/D  =  12,  r/D  =  0.5,  as  shown  in  Fig.  40. 

The  shapes  and  positions  of  each  interface  were  also  recorded 
at  intervals  of  16  ms  after  passing  through  the  trigger  position. 
Results  obtained  by  using  this  procedure  are  shown  in  Fig.  41. 

The  averaged  shapes  of  the  interfaces,  at  t  =  0  are  similar  to 
the  potential/irrotational  flow  interfaces  found  in  cold  jets. 

The  probability  distribution  of  interface  positions,  at 
r/D  =  0.5;  x  =  16  ms,  has  a  totalled  measured  range  of  positions 
of  plus  or  minus  12%  about  the  mean.  It  is  expected  that  this 
variation,  coupled  with  variation  in  interface  structures,  will 
contribute  to  an  appreciable  phase  scrambling  effect  on  multi¬ 
probe  measurements .  The  double  ionization  probe  measurements 
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which  have  been  described  will  provide  a  comparison  with  the 
film  analysis  data. 
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It  has  been  found  that  coherent  structures  are  at  least 
as  important  in  jet  flames  as  they  are  in  cold  jets.  The 
experiments  now  in  progress  will  provide  quantitative  information 
and  descriptions  of  these  structures  in  jet  flames.  Furthermore, 
the  experiments  will  be  extended  downstream  to  investigate  the 
full  length  of  the  flame.  An  important  question  which  still 
remains  unanswered  is  the  extent  to  which  these  flames  can  be 
considered  to  be  turbulent.  At  x  =  8  D  the  flame  is  undoubtedly 
transitional,  containing  periodic,  essentially  axisymmetric, 
vortex  rings .  Beyond  x  =  14  D,  the  flame  has  all  the  appearances 
of  turbulent  three-dimensional  flow,  although  coherent  structures 
are  still  present. 
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TABLE  1 


Jet  Flame  I  Parameters 


Primary  (diluting)  air  flow  rate 
Primary  propane  flow  rate 
Equivalence  ratio 
Uj,  primary  jet  velocity 
Ug,  secondary  flow  velocity 

Cold  Jet  1  Parameters 


Primary  air  flow  rate 
Primary  propane  flow  rate 
Uj 
US 

<Sj,  Nozzle  boundary  layer  (shear)  thickness 

Nozzle  Geometry 

D,  nozzle  internal  diameter 

Secondary  flow  diameter  (square  cross-section) 
Nozzle  lip  thickness 

Physical  Gas  Properties  at  2o°c  and  1  at ms . 

Propane  (C3Hq)  density 
Propane  viscosity 
Air  density 
Air  viscosity 


120  1/min 
44  1/min 

10.4 

6.4  m/s 
0.74  m/s 


164  1/min 
0 

6.4  m/s 
0.74  m/s 
2 .36  mm 


25.4  mm 
400  mm 
1  mm 


y.3  kg/m3 
80.29  u  poise 
1.205  kg/m3 
183  y  poise 


Mean  velocity 
profiles 


Cartesian 
coordinates 
based  on 
nozzle  lip 


U  +  u 


Velocity 

components 


Cylindrical 

coordinates 


Fig.  1  Notation  for  axisymmetric  jet  flame 


developed  region 


ri9.  2  The  physical  structure  in  ths  nixing  layer 
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Fig.  3  Apparatus  for  the  investigation  of  axisymmetric  jet  flame 


fig.  4  Laser  Doppler  Anemometer  System 


Scattering  Particle 
Direction 


Doppler  signal  -  single  scattering  particle 


counting  sytt< 


Fig.  7  Designs  of  fine  wire  thermocouples 


System  for  measuremment  of  thermocouple  time  constant 
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DIMENSIONS  IN  mm 
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Fig.  12  Ionization  probe  circuit 


Fig.  13  Colour  Schlieren  system 


Fig.  14  Schematic  view  of  transitional  flame 


IONIZATION  PROBE 


IONIZATION  PROBE 


using  Two  Ionization  Probes 


Phase  Scrambling  Effects  on  Signal  Recovery  in  Conditional  Sampling 


Fig.  19a  Photograph  of  Flame  I 
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Fig.  24  Nozzle  boundary  layer  thickness  versus  Reynolds  number 


Fig.  26  Centreline  distributions  of  mean  velocity  and  turbulence  intensity  in  Cold  Jet  I  with  and 
without  a  secondary  flow 


Fig.  27  Amplitude  spectra  of  u,  Cold  Jet 


Fig.  28  Amplitude  spectra  of  u.  Cold  Jet  I,  r  =  0#  2  <  x/D 


Fig.  32  Amplitude  spectra  of  u.  Cold  Jet  I,  x/D 
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Fig.  34  Comparison  of  burning  and  nonburning  jets  at  x/d  =  4 
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Fig.  36  Amplitude  spectrum  from  ionization  probe,  x/D  =  8;  r/D 
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Fig.  37  Comparison  of  fundamental  instability  frequency  in  Cold  Jet  I  with  theory  of 
Michalke 
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Fig.  40  Schematic  view  of  Jet  flame  showing  visual  appearance  of  high-temperature 
interfaces  and  the  method  for  ensemhie  averaging  from  films 


